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Abstract 

We present a comprehensive overview of different types of parametric interactions in nonlinear optics wliicli 
are associated witfi simultaneous phase-matching of several optical processes in quadratic nonlinear media, 
the so-called multistep parametric interactions. We discuss a number of possibilities of double and multiple 
phase-matching in engineered structures with the sign-varying second-order nonlinear susceptibility, including 
(i) uniform and non-uniform quasi-phase-matched (QPM) periodic optical superlattices, (ii) phase-reversed and 
periodically chirped QPM structures, and (iii) uniform QPM structures in non-coUinear geometry, including re- 
cently fabricated two-dimensional nonlinear quadratic photonic crystals. We also summarize the most important 
experimental results on the multi-frequency generation due to multistep parametric processes, and overview the 
physics and basic properties of multi-color optical parametric solitons generated by these parametric interactions. 
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§ 1. Introduction 

Energy transfer between different modes and phase-matching relations are the fundamental concepts in nonlinear 
optics. Unlike nonparametric nonlinear processes such as self-action and self-focusing of light in a nonlinear Kerr- 
like medium, parametric processes involve several waves at different frequencies and they require special relations 
between the wave numbers and wave group velocities to be satisfied, the so-called phase-matching conditions. 

Parametric coupling between waves occurs naturally in nonlinear materials without the inversion symmetry, when 
the lowest-order nonlinear effects are presented by quadratic nonlinearities, often called x^'^'^ nonlinearities because 
they are associated with the second-order contribution (^^ x'^^^^) to the nonlinear polarization of a medium. 
Conventionally, the phase-matching conditions for most parametric processes in optics are implemented either by 
using anisotropic crystals (the so-called perfect phase-matching), or they occur in fabricated structures with a 
periodically reversed sign of the quadratic susceptibility (the so-called quasi-phase-matching or QPM). The QPM 
technique is one of the leading technologies these days, and it employs the spatial scales (^^ 1 -f- which are 

compatible with the operational wavelengths of optical communication systems. 

Nonlinear effects produced by quadratic intensity-dependent response of a transparent dielectric medium are 
usually associated with parametric frequency conversion such as the second harmonic generation (SHG). The SHG 
process is one of the most well-studied parametric interactions which may occur in a quadratic nonlinear medium. 
Moreover, recent theoretical and experimental results demonstrate that quadratic nonlinearities can also produce 
many of the effects attributed to nonresonant Kerr nonlinearities via cascading of several second-order parametric 
processes. Such second-order cascading effects can simulate third-order processes and, in particular, those associ- 
ated with the intensity-dependent change of the medium refractive index Stegeman, Hagan, and Torner [1996]. 
Importantly, the effective (or induced) cubic nonlinearity resulting from a cascaded SHG process in a quadratic 
medium can be of the several orders of magnitude higher than that usually measured in centrosymmetric Kerr-like 
nonlinear media, and it is practically instantaneous. 

The simplest type of the phase-matched parametric interaction is based on the simultaneous action of two second- 
order parametric sub-processes that belong to a single second-order interaction. For example, the so-called two-step 
cascading associated with type I SHG includes the generation of the second harmonic (SH), oj + oj ^ 2uj, followed 
by the reconstruction of the fundamental wave through the down-conversion frequency mixing process, 2uj~uj ^ uj. 
These two sub-processes depend only on a single phase-matching parameter Ak. In particular, for the nonlinear 
X*^^) media with a periodic modulation of the quadratic nonlinearity, for the QPM periodic structures, we have 
Afc ~ k2 — 2fci + Gm, where ki = k{uj), k2 = k{2iu) and Gm is a reciprocal vector of the periodic structure, 
Gm = 277111/ A, where A is the lattice spacing and m is integer. For a homogeneous bulk x^^-* medium, we have 

Gm = 0. 

Multistep parametric interactions and multistep cascading presents a special type of the second-order parametric 
processes that involve several different second-order nonlinear interactions; they are characterized by at least two 
different phase-matching parameters. For example, two parent processes of the so-called third-harmonic cascading 
are: (i) second-harmonic generation, uj + lu = 2uj, and (ii) sum-frequency mixing, to -\~ 2uj — 3a;. Here, we may 
distinguish five harmonic sub-processes, and the multistep interaction results in their simultaneous action. 

Different types of the multistep parametric processes include third-harmonic cascaded generation, two-color para- 
metric interaction, fourth-harmonic cascading, difference-frequency generation, etc. Various applications of the 
multistep parametric processes have been mentioned in the literature. In particular, the multistep parametric in- 
teraction can support multi-color solitary waves, it usually leads to larger accumulated nonlinear phase shifts, in 
comparison with the simple cascading, it can be effectively employed for the simultaneous generation of higher-order 
harmonics in a single quadratic crystal, and also for the generation of a cross-polarized wave and frequency shifting 
in fiber optics gratings. Generally, simultaneous phase-matching of several parametric processes cannot be achieved 
by the traditional methods such as those based on the optical birefringence effect. However, the situation becomes 
different for the media with a periodic change of the sign of the quadratic nonlinearity, as occurs in the QPM 
structures or two-dimensional nonlinear photonic crystals. 

In this review paper, we describe the basic principles of simultaneous phase-matching of two (or more) parametric 
processes in different types of one- and two-dimensional nonlinear quadratic optical lattices. We divide different 
types of phase-matched parametric processes studied in nonlinear optics into two major classes, as shown in Fig. 1, 
and discuss different types of parametric interactions associated with simultaneous phase-matching of several optical 
processes in quadratic (or x^^') nonlinear media, the so-called multistep parametric interactions. In particular, 
we overview the basic principles of double and multiple phase-matching in engineered structures with the sign- 
varying second-order nonlinear susceptibility, including different types of QPM optical superlattices, non-coUinear 
geometry, and two-dimensional nonlinear quadratic photonic crystals (which can be considered as two-dimensional 
QPM lattices). We also summarize the most important experimental results on the multi- frequency generation 
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Figure 1: Different types of parametric processes in nonlinear optics, and the specific topics covered by this review 
paper. SHG: second-harmonic generation; SFG and DFG: sum- and difference-frequency generation. THG: third- 
harmonic generation; FWM: four-wave mixing; DFWM: degenerate FWM; OPO: optical parametric oscillator. 



due to multistep parametric processes, and overview the physics and basic properties of multi-color optical solitons 
generated by these parametric interactions. 

§ 2. Single-phase-matched processes 

One of the simplest and first studied parametric process in nonlinear optics is the second-harmonic generation 
(SHG). The SHG process is a special case of a more general three- wave mixing process which occurs in a dielectric 
medium with a quadratic intensity-dependent response. The three-wave mixing and SHG processes require only 
one phase-matching condition to be satisfied and, therefore, they both can be classified as single phase-matched 
processes. 

In this section, we discuss briefly these single phase-matched processes, and consider parametric interaction 
between three continuous-wave (CW) waves with the electric fields Ej = ■^[Aj exp(— ikj • r + iojjt) + c.c], where 
j = 1,2,3, with the three frequencies satisfying the energy-conservation condition, lui + UJ2 = W3. We assume that 
the phase-matching condition is nearly satisfied, with a small mismatch Afc among the three wave vectors; i.e., 
Afc = k^lius) — ki{uji) — fc2(^2)- In general, the three waves do not propagate along the same direction, and the 
beams may walk off from each other as they propagate inside the crystal. If all three wave vectors point along the 
same direction (as, e.g., in the case of the QPM materials), the waves have the same phase velocity and exhibit no 
walk-off. 

The theory of x^^^-mediated three-wave mixing is available in several books devoted to nonlinear optics (Shcn 
[1984]; Butcher and Cotter [1992]; Boyd [1992]). The starting point is the Maxwell wave equation written as 

1 d^E 1 d^F , , 

VxVxE+-,^ = --,^, (2.1) 

where eg is the vacuum permittivity and c is the speed of light in a vacuum. The induced polarization is written in 
in the frequency domain as 

P(r, Lo) = eoX^'^E + eox^'^EE + • • • , (2.2) 
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where a tilde denotes the Fourier transform. Using the slowly varying envelope approximation and neglecting the 
walk-off, one can derive the following set of three coupled equations describing the parametric interaction of three 
waves under type II phase matching: 

0, (2.3) 
0, (2.4) 
0, (2.5) 

where F = d^J /c^, and d^g is a convolution of the second-order susceptibility tensor x^'^' and the polarization unit 

vectors of the three fields, d^g = I < e3x'-^-'eie2 >. 

In the case of type I SHG, only a single beam at the pump frequency tui is incident on the nonlinear crystal, and 
a new optical field at the frequency 2loi is generated during the SHG process. We can adapt Eqs. (2.3)-(2.5) to 
this case with minor modifications. More specifically, we set tu^ = 2lui and Ai = A2. The first two equations then 
become identical, and one of them can be dropped. The type I SHG process is thus governed by the following set 
of of two coupled equations: 

0, (2.6) 
0, (2.7) 

where a = {uji / nic)d^^^ is the nonlinearity parameter and Afc = ^3 — 2ki is the phase-mismatch parameter. 

Both the three- wave mixing and SHG processes present an example of a singlc-phase-matched parametric process 
because it is controlled by a single phase- matching parameter Afc. This kind of parametric processes can be described 
as a two-step cascading interaction which includes: (i) the generation of the second-harmonic (SH) wave, u+uj = 2a;, 
followed by (ii) the reconstruction of the fundamental wave through the down-conversion frequency mixing process, 
i.e. 2uj ~ Lo = uj. Respectively, the first sub-process is responsible for the generation of the SH field, with the 
most efficient conversion observed at Afc = 0, while the second sub-process, also called cascading, can be associated 
with an effective intensity-dependent change of the phase of the fundamental harmonic (^ d^/Afc), which is similar 
to that of the cubic nonlinearity (DeSalvo, Hagan, Sheik-Bahae, Stegeman, Vanstryland, and Vanherzeele [1992]; 
Stegcman, Hagan, and Torncr [1996]; Assanto, Stegeman, Sheik-Bahae, and Vanstryland [1995]). This latter 
effect is responsible for the generation of the so-called quadratic solitons, two- wave parametric soliton composed of 
the mutually coupled fundamental and second-harmonic components (Sukhorukov [1988]; Torner [1998]; Kivshar 
[1997]; Etrich, Lcderer, Malomed, Pcschel, and Peschel [2000]; Torruellas, Kivshar, and Stegeman [2001]; Boardman 
and Sukhorukov [2001]; Buryak, Di Trapani, Skryabin, and Trillo [2002], and references therein). 

Multistep parametric interactions and multistcp cascading effects discussed in this review paper are presented 
by different types of phase-matched parametric interactions in quadratic (or x'^') nonlinear media which involve 
several different parametrically interacting waves, e.g. as in the case of the frequency mixing and sum-frequency 
generation. However, all such interactions can also be associated with the two major physical mechanisms of the 
wave interaction discussed for the SH process above: (i) parametric energy transfer between waves determined by 
the phase-mismatch between the wave vectors of the interacting waves, and (ii) phase changes due to this parametric 
interaction. Below, we discuss these interactions for a number of physically important examples. 

§ 3. Multistep phase-matched interactions 



dz 

tfca^-c^iFAaAte-''^^- = 
dz 

,k,'^-ulTA,A2e^'^^ = 



dz 

ik, ^ _ 2k,cTAy^''' 
dz 



Table 1: Examples of multistep parametric interactions involving SHG. 

No. Multistep parametric pro- Cascading x^^^ Equivalent high- 
ccss steps order parametric 

process 

1 Type I third-harmonic mul- lj + cj — 2lj; lu -\- u) ^ u) = 'iuj 
tistep process 2uj — Sco 

2 Type II third-harmonic mul- uj + co = 2uj ; uj± + uj + uj + uj± = 3a; 
tistep process 2aj = 3a; 
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No. Multistep parametric pro- Cascading x'^^ Equivalent high- 
cess steps order parametric 

process 

3 3:1 frequency conversion and 'iuj 2uj + u \ u = 
division 2lu — uj 

4 Fourth-harmonic multistep uj -\- u) = 2uj \ 2uli -\- uj-\-uj-\-uj-\-uj = Auj 
process 2uj — 4ui 

5 Type I & Type II two-color uj + lu — 2uj \ 2uj — uj ^ uj — u = 
multistep process = ojx 

6 Type I & Type I two-color \ u) = 2u} \ iox_ + 
multistep process u)x_ = 2ct; 

7 wavelength conversion -f = 2aj ; 2a; — u ^- — iDa = 

8 Self-doubling OPO ujp uji + ujs ; <jJs + 

9 Self-sum- frequency genera- ujp u)i-\~ujs ; ujs-\- 
tion OPO LOp = tJsPG 

10 Internally pumped OPO + ^p/2 = \ 



In the table above, the symbol lj±_ stands for a wave polarized in the plane perpendicular to that of the wave with the main 
carrier frequency uj. 

In this part, we consider the nonlinear parametric interactions that involve several processes, such that each of 
the processes is described by an independent phase-matching parameter. In the early days of nonlinear optics, the 
motivation to study this kind of parametric interactions was to explore various possibilities for the simultaneous 
generation of several harmonics in a single nonlinear crystal (see, e.g., Akhmanov and Khokhlov [1964, 1972]) as 
well as to use the cascading of several parametric processes for measuring higher-order susceptibilities in nonlinear 
optical crystals (see, e.g., Yablonovitch, Flytzanis, and Bloembergen [1972]; Akhmanov, Dubovik, Saltiel, Tomov, 
and Tunkin [1974]; Akhmanov [1977]; Kildal and Iseler [1979]; Bloembergen [1982]). More recently, these 
processes were proved to be efficient for the higher-order harmonic generation, for building reliable standards for 
the third-order nonlinear susceptibility measurements (see, e.g., Bosshard, Gubler, Kaatz, Mazerant, and Meier 
[2000] and references therein), and also for generating multi-color optical solitons. Additionally, it is expected that 
the multistep parametric processes and multistep cascading will find their applications in optical communication 
devices, for wavelength shifting and all-optical switching (see discussions below). Another class of applications of 
the multistep phase-matched parametric processes is the construction of optical parametric oscillators (OPOs) and 
optical parametric amplifiers (OPAs) with complimentary phase-matched processes in the same nonlinear crystal 
where the main phase-matched parametric process occurs. In this way, the OPOs and OPAs may possess additional 
coherent tunable outputs (see Sec. 3.5 below). On the basis of the third-harmonic multistep cascading process 
{w + UJ = 2u>] UJ -\- 2uj = 'iiij } , real advances have been made in the development of the nonlinear optical systems 
for division by three. The multistep parametric interactions governed by several phase-mismatching parameters 
can also occur in centrosymmetric nonlinear media (see, e.g., Akhmanov, Martynov, Saltiel, and Tunkin [1975]; 
Reintjes [1984]; Astinov, Kubarych, Milne, and Miller [2000]; Crespo, Mendonca, and Dos Santos [2000]; Misoguti, 
Backus, Durfee, Bartels, Murnane, and Kapteyn [2001]). 

The first proposal for simultaneous phase-matching of two parametric processes can be found in the pioneering 
book of Akhmanov and Khokhlov (Akhmanov and Khokhlov [1964, 1972]) who derived and investigated the 
condition for the THG process in a single crystal with the quadratic nonlinearity through the combined action of the 
SHG and SFG parametric processes. For the efficient frequency conversion, both the parametric processes, i.e. SHG 
and SFG, should be phase matched simultaneously. Some earlier experimental attempts to achieve simultaneously 
two phase-matching conditions for the SHG and SFG processes were not very successful (Sukhorukov and Tomov 
[1970], Orlov, Sukhorukov, and Tomov [1972]). Recent development of novel techniques for efficient phase matching, 
including the QPM technique, make many of such multistep processes readily possible. 

The multistep parametric processes investigated so far can be divided into several groups, as shown in Fig. 1 
and Table 1. These processes include: third-harmonic multistep cascading; fourth- harmonic multistep cascading; 
two-color multistep cascading; FWM and DFWM in x*^^' media , OPO and OPA multistep cascading. Wc will 
review these groups separately. 
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Figure 2: Schematic of the THG multistep cascading. 



Table 1 does not provide a complete list of all possible types of the multistep parametric processes. Here, 
we pick up only several examples for each type of the multistep process. As follows from Table 1, some of the 
parametric processes simulate some known higher-order processes but occur through several steps. However, some 
other multistep parametric processes have no such simpler analogues. Below in the review paper, we discuss some 
of those processes in more details. For example, the THG multistep cascading is considered in Sec. 3.1, while 
the two-color multistep cascading is discussed in Sec. 3.3. The last three lines in Table 1 are the examples of the 
multistep cascading processes in optical parametric oscillators and amplifiers; and they will be discussed in Sec. 3.5. 

To illustrate the physics responsible for the use of the terminology "multistep interaction" or "multistep cascad- 
ing", we just point out the simultaneous action of SHG and SFG (line 1 in Table 1); this three- wave interaction 
involves five simpler parametric sub-processes in a quadratic medium. In order to describe, for example, the non- 
linear phase shift of the fundamental wave accumulated in this interaction, we should consider the following chain 
of parametric interactions: SHG {uj + ui = 2tj), SFG {cu + 2u! = Sw), DFM (3cj — uj = 2w), and, finally, another 
DFM {2uj-uj = uj). 

3.1. THIRD-HARMONIC MULTISTEP PROCESSES 

The multistep parametric interaction involving the THG process is one of the most extensively studied multistep 
cascading schemes (see Fig. 2). The two simpler parametric interactions are the SHG process, uj + oj = 2ijj, and 
the SFG process, uj + 2lu ^ 3uj. Each of these sub-processes is characterized by an independent phase matching 
parameter, namely A/cshg = k2uj — 2fciu and A/csfg = k^^j — — few 

This kind of the multistep cascading appears in many schemes of parametric interactions in nonlinear optics, 
including: (i) efficient generation of the third harmonic in a single quadratic crystal; (ii) measurement of unknown 
X^^^ tensor components using known x^^^ components of the crystal as a reference; (iii) accumulation of a large 
nonlinear phase shift by the fundamental wave; (iv) propagation of multi-color solitons; (v) frequency division; (vi) 
generation of entangled and squeezed photon states, etc. 

3.1.1. Efficient generation of a third- harmonic wave 

First we consider the process of efficient generation of a third-harmonic wave. In the approximation of plane waves, 
this process is described by the following system of coupled equations for the slowly varying amplitudes Ai, A2, 
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and A3 of the fundamental, second-, and third-harmonic waves, respectively: 

-iAfcsFG^ 



dAi 
dz 

dA2 
dz 

dAs 
dz 



-iaiA2^Ie-*'^'=^»°" - ia3A3A;e' 



-ia5A2Aie'^''<'''''' - ijAle'^''™'^'' , 



(3.1) 



fci + Gq, (71,2 = (27r/Aini,2)dcff,i and aj 



where A/cshg = — 2/ci + Gp and A/csfg = ^3 — 
{ujj-2/i-ui){2TT/ Xinj-2)dcsji (where j — 3,4,5). Here, the parameters c?off,i and deft, 11 are the effective quadratic 
nonlinearities corresponding to the two steps of the multistep parametric process; its values depend on the crystal 
orientation (see, e.g., Dmitricv et al. [1999]) and the method of phase matching (PM). The parameter 7 is found as 
7 = (37r/4Ai7i3)Xoff (for calculation of in crystals see Yang and Xie [1995]). The complementary wave vectors 
Gp and Gq are two vectors of the QPM structure that can be used for achieving double phase matching. A solution 
of this system (see, e.g., Kim and Yoon [2002]; Qin, Zhu, Zhang, and Ming [2003]) gives maximum for THG in sev- 
eral different situations, when (i) AfcsHG — >■ 0; (ii) A/csfg — *■ 0; (iii) A/cthg = A^shg + A/csfg = — 3/ci — *■ 0; 
and (iv) simultaneously AfcsHG — > and AfcspG — ^ 0. The latter condition, for which both SHG and SFG 
steps should be simultaneously phase matched, corresponds to the highest efficiency for THG. The intensity of the 
third-harmonic (TH) wave is proportional to the fourth power of the crystal length. 



\A3{L)\^^\alal\A,\'^L\ 
This expression should be compared with the analogous result for the ccntrosymmetric media: 

1^3(i)|' = 7'l^i|'i'- 



(3.2) 



(3.3) 



The advantage of the single-crystal phase-matched cascaded THG becomes clear if we note that, in average, {(J2<JzL)^ 
is in 10"* — 10^ times larger than 7^, even for the sample length as small as 1 mm. 

Table 2: Experimental results on cascaded THG processes 



Nonlinear 




Phase-matched steps 


Phase matching 


L 


Regime 


V,% 


Refs. 


crystal 


[^m] 






method 


[cm] 








LiTaOa 


1.44 


AfcsHG 


^ 0; AfcsFG ~ 


QPOS 


1.5 


Pulsed 
(8 ns) 

Pulsed 
(8 ns) 

quasi- cw 
(30 ns) 


27% 


a 


LiTaOs 


1.570 


AfcsHG 


^ 0; AfcsFG ~ 


QPOS 


0.8 


23% 


b, c 


LiTaOa 


1.342 


AfcsHG 


^ 0; AfcsFG ~ 


SHG-QPM Ist-ord. 
SFG-QPM 3rd-ord. 


1.8 


19.3 % 


d 


LiTaOa 


1.342 


A/CSHG 


^ 0; AfcsFG ~ 


SHG-QPM Ist-ord. 
SFG-QPM 3rd-ord. 


1.2 


Pulsed 
(90 ns) 


10.2% 


e 


/3-BBO 


1.055 


AfcxHG 
A/csFG 


= -AfcsHG / 


BPM 


0.3 


Pulsed 
(350 fs) 


6% 


f,g 


LiTaOa 


1.442 


AksHG 


^ 0; AfcsFG ~ 


QPOS 


0.6 


Pulsed 
(10 ns) 


5.8% 


h 


KTP 


1.8 


AfcsHG 


/ 0; AfcsFG ~ 


BPM 




Pulsed 
(35 ps) 


5% 


i 


LiTaOa 


1.064 


AfcsHG 


^ 0; AfcsFG ~ 


PR-QPM 


1.2 


quasi- cw 
(150 ns) 


2.8% 


j 


KTP 


1.618 


AfcTHG 
AfcsFG 


= -AfcsHG / 


BPM 


0.11 


Pulsed 
(22 ps) 


2.4% 


k, 1 


SBN 


1.728 


AfcsHG 


^ 0; AfcsFG ~ 


QPOS 


0.75 


Pulsed 
(15 ps) 


1.6% 


m 


d-LAP 


1.055 


AfcTHG 
AfcsFG 


= -AfcsHG / 


BPM 


0.1 


Pulsed 

(350 fs) 


1.2% 
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Nonlinear 



crystal 



Phase-matched steps 



[^^m] 



Phase matching 
method 



Regime r;,% 



Refs. 



/3-BBO 
KTP 

waveguide 
KTP 

LiNbOs 
waveguide 

LiNbOs 
KTP 

waveguide 
/3-BBO 

LiNbOs 
2D-NPC 

LiNbOa 
KTP 

YiLiNbOa 



1.05 AfcxHG ~ 0; 

AfcsFG = — AfcsHG / 

L234 AksHG ~ 0; AksFc ~ 

1.32 AfcsHG / 0; AfcsFG ~ 

1.619 AfcsHG ~ 0; AfcsPG ~ 

1.534 AfcsHG / 0; AfcsPG ~ 

1.65 AfcsHG ~ 0; AfcsFG ~ 

1.053 AfeTHG « 0; 

AfesFG = -AfcsHG / 

1.536 AfcsHG ~ 0; AfcsFG ~ 

3.561 AfcsHG ~ 0; AfcsPG ~ 

1.55 AfcsHG « 0; AfcsFG « 

1.064 AfcsHG ~ 0; AfcsFG ~ 



BPM 
QPM 
BPM 

SHG- QPM Ist-ord. 
SFG- QPM 3rd-ord. 

QPM 

SHG- QPM Ist-ord. 
SFG- QPM 3rd-ord. 

BPM 
2D-QPM 
QPM 
QPOS 

SHG-QPM 9th-ord. 
SFG-QPM 33rd-ord. 



0.72 Pulsed 
(5 ps) 

0.26 Pulsed 
(9 ns) 

0.47 pulsed 
(200 fs) 

pulsed 
(7ps) 

0.6 Pulsed 
(9 ns) 

0.35 pulsed 
(6 ps) 

0.7 Pulsed 
(45 ps) 

1 pulsed 
(5 ns) 

2 CW 

1 CW 

0.5 pulsed 
(100 ns) 



0.8% 

0.4%2 

0.17% 

0.055% 

0.016%^ 

0.011% 

0.007% 

0.01%^ 

10-*% 

3.10-^% 
[W.cm-2] 

10-5% 



p,q 
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Mu, Gu, Makarov, Ding, Wang, Wei, and Liu [2000] 
" Ding, Mu, and Gu [2000] 

Baldi, Trevino-Palacios, Stegeman, Demicheli, Ostrowsky, Delacourt, and Papuchon [1995] 
" Gu, Korotkov, Ding, Kang, and Khurgin [1998] 
' Sundheimer, Villeneuve, Stegeman, and Bierlein [1994b] 
" Tomov, Van Wonterghem, and Rentzepis [1992] 
" Broderick, Bratfalean, Monro, Richardson, and de Sterke [2002] 
™ Pfister, Wells, HoUberg, Zink, Van Baak, Levenson, and Bosenberg [1997] 
^ Fradkin-Kashi, Arie, Urenski, and Rosenman [2002] 



2 backward THG 

3 backward THG 
*Type II cascaded THG 
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Figure 3: (above) Schematic of the experimental setup, (below) Average powers of the second harmonic (right 
scale) and third harmonic (left scale) vs. temperature. Average power of the fundamental wave is 4.8 mW (Zhang, 
Wei, Zhu, Wang, Zhu, and Ming [2001]). 



^ Volkov, Laptev, Morozov, Naumova, and Chirkin [1998] 
Abbreviations: 

QPOS - quasi periodical optical superlattices; BPM - birefringence phase matching; PR-QPM - phase reversed QPM 
structure (see Chou, Paramcswaran, Fejer, and Brener [1999]); 2D-NPC - two-dimensional photonic crystals; 2D-QPM - 
two-dimensional QPM structure. 



Introducing normalized efficiency (measured in units of W ^cm ^) for the first and second steps in separate 
crystals as 770,1 and 770,2, respectively, we obtain 

= lvo,im,2P?L^ (3.4) 

The results for THG in a single quadratic crystal under the condition of double phase matching were reported for 
the efficiency exceeding 20%, as shown in Table 2. Figure 3 shows the phase matching curves in the experiment 
(Zhang, Wei, Zhu, Wang, Zhu, and Ming [2001]) where 27% THG efficiency has been achieved. The two PM curves 
are not perfectly overlapped. We may expect that with an improvement of the superlattice structure the achieved 
efficiencies will be higher. Numerical solution of the system (3.1) shows that, in general, the third-harmonic output 
has an oscillating behavior as a function of the length or input power, and generally it does not reach the efficiency 
of 100%. However, as shown in Egorov and Sukhorukov [1998]; Chirkin, Volkov, Laptev, and Morozov [2000] and 
Zhang, Zhu, Yang, Qin, Zhu, Chen, Liu, and Ming [2000], the total conversion of the fundamental wave into the 
third-harmonic wave is possible if the ratio 02 1 '^^ is optimized. 

In other cases, where only one phase-matched condition is satisfied, the THG conversion efficiency is not so large, 
and it should be compared to that of the direct THG process in a cubic medium. The efficiency of a cascaded THG 
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process is inversely proportional to the square of the wave vector mismatch of the unmatched process. It is important 
to note that in such cases, even when only one of the phase-matched parameters AfcsHG or AfcsFG is close to zero 
(such conditions can be easily realized in the birefringcnt phase matching in bulk quadratic media or in an uniform 
QPM structure) the THG process has the behavior of the phase matched process with the characteristic dependence 
[sin(x)/a;]^ for the TH intensity vs. tuning and a cubic dependence on the input intensity. The dependence of the 
sample thickness is quadratic but not periodic function as for the totally non phase matched processes. 

Situation (iii), that is, tuning where A/cshg + A/csfg = 0, also corresponds to the phase-matching condition for 
direct THG where the fundamental wave is converted directly into a third-harmonic wave. In such a case, both the 
cascade process and direct process = 3fci contribute to the third harmonic wave. However, a relative contribution 
of the two processes could be different. In some cases, the direct THG process is stronger (see, e.g., Feve, Boulanger, 
and Guillien [2000]), in other cases the cascading THG process is dominant (see, e.g.. Banks, Fcit, and Perry 
[2002]; Bosshard, Gubler, Kaatz, Mazcrant, and Meier [2000]). 

We wonder why the whole process is phase matched when both steps are mismatched? The reason is that in 
the case (iii) we have the situation similar to that of the quasi-phase- matching effect (see, e.g., Reintjcs [1984]; 
Banks, Fcit, and Perry [1999]; Durfec, Misoguti, Backus, Kapteyn, and Murnane [2002]). Indeed, if the first step 
is mismatched, inside a nonlinear medium it generates a periodically modulated polarization at the frequency 2u; 
(P(2aj)) and the modulation period is exactly a mismatch of the second step. Thus, in the regions where the phase 
of the polarization P(2omcga) is reversed, we have the minimum generation of second- and third- harmonic waves, 
and the generated TH components interfere constructively with the propagating third-harmonic wave along the 
length of the crystal. This leads to a quadratic dependence of the THG efficiency on the crystal length. 

For the cascaded THG processes, the optimal focusing is an important issue when the goal is the maximum 
conversion efficiency. If only one of the steps is phase matched, the optimal focusing is in the input face when SFG 
is phase matched, or at the output face when SHG is phase matched (Rostovtscva, Sukhorukov, Tunkin, and Salticl 
[1977]; Rostovtscva, Saltiel, Sukhorukov, and Timkin [1980]). If both the steps are phase matched, the optimum 
focusing position is in the center of the nonlinear medium (Ivanov, Koynov, and Salticl [2002]). 

In Table 2, wc summarize the experimental results obtained for the efficient single-crystal THG processes. Con- 
ditions at which the third-harmonic wave included in double phase-matched interaction can be transformed with 
100% efficiency into the 2uj wave or uj wave were found by Komissarova and Sukhorukov [1993]; Komissarova, 
Sukhorukov, and Tcrcshkov [1997]; and Egorov and Sukhorukov [1998]. Volkov and Chirkin [1998] show that 
the same parametric interaction can be used for 100% conversion from 2ll} wave into Suj wave. Considering the 
situations with nonzero SH and TH boundary conditions it has to be taken into account that, as shown by Alekseev 
and Ponomarcv [2002], the spatial evolution of three light waves participating simultaneously in SHG and SFM 
under the conditions of QPM double phase matching becomes chaotic at large propagation distances for many 
values of the complex input wave amplitudes. Thus, the possibility of transition to chaos exists in the application of 
an additional pump at frequency 2uj in order to increase the efficiency of THG [Egorov and Sukhorukov [1998]]. As 
shown in Longhi [2001b] for the multistep cascading process {cu + uj = 2uj; w -I- 2w = So; } in a cavity, the formation 
of spatial patterns is possible. Another application of the multistep cascading interaction and, in particular, the 
THG multistep cascading is the generation of the entangled and squeezed quantum states (see Sec. 3.5 below). The 
conditions for the simultaneous phase matching of both SHG and SFG processes were considered by Pfister, Wells, 
HoUberg, Zink, Van Baak, Lcvenson, and Bosenberg [1997]; Grechin and Dmitriev [2001a], for uniform QPM 
structures; by Zhu and Ming [1999], for quasi-periodic optical superlattices; by Fradkin-Kashi and Arie [1999] and 
Fradkin-Kashi, Aric, Urcnski, and Roscnman [2002], for the generalized Fibonacci structures; and by Saltiel and 
Kivshar [2000a], for the two-dimensional nonlinear photonic crystals. This topic will be discussed below in Sec. 4. 

3.1.2. Nonlinear phase shift in multistep cascading 

The multistep parametric process, that combines two phase-matched interaction SHG (uj + lu ~ 2uj) and SFG 
{2lu + lj = 3cj), is described by the system of equations (3.1), and it can be used for all-optical processing and 
formation of parametric solitons. Indeed, the nonlinear phase shift (NPS) accumulated by the fundamental wave 
in the multistep cascading process is in several times larger than that in the standard cascading interaction that 
involves only one phase matched process (Koynov and Saltiel [1998]). 

To illustrate this result, we consider the fundamental wave with the frequency uj entering a second-order nonlinear 
media under appropriate phase-matching conditions. As the first step, the wave with frequency 2uj is generated via 
the type I SHG process, and then, as the second step, the third-harmonic wave is generated via the SFG process 
{2iij + uj = 3u!). Both the processes, SHG and SFG, are assumed to be nearly phase matched. The generated second- 
and third- harmonic waves arc down-converted to the fundamental wave to via the processes (2a; — to); (3a; — 2lu), 
and (3w — w, 2a; — w), contributing all to the nonlinear phase shift that the fundamental wave collects. As follows 
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Figure 4: Schematic of the possible channels for the phase modulation of the fundamental wave (Koynov and Saltiel 
[1998]). 
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Figure 5: Nonlinear phase shift and depletion of the fundamental wave as a function of its normalized input 
amplitude: solid line - multistep cascading, dash line - type I SHG case (Koynov and Saltiel [1998]). 



from numerical calculations, the total NPS is a result of the simultaneous action of two-, three-, and four-step 
X*-^^ cascading, and it can exceed the value of tt for relatively low input intensities. The possible channels for the 
phase modulation and NPS of the fundamental wave are shown in Fig. 4. The interpretation of the analytical 
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Figure 6: Wavelength conversion of optical communication channels in a periodically-poled nonlinear crystal using 
SHG/DFG multistep cascading. 



result obtained in the fixed- intensity approximation (Koynov and Saltiel [1998]) shows that the effective cascaded 
fifth-order and higher-order nonlinearities are involved into the accumulation of a total NFS, and the signs of 
the contributions from different processes can be controlled by a small change of the phase matching conditions. 
Schematically, the role of the multistep cascading with three and four steps can be interpreted as equivalent of a 
contribution from the higher-order nonlinear corrections 714 and rig to the refractive index, and they can be linked 
to the cascaded x^^-' and x^^-* processes, respectively. Indeed, for a relatively low intensity of the fundamental wave, 
the refractive index can be written in the form of expansion, 

n{E) = no + n-zE^ + mE* + naE^ + .... (3.5) 

The advantage of the multistep cascading for accumulating large nonlinear phase shift over the conventional 
two-step cascading is illustrated in Fig. 5 obtained by numerical integration of the system (3.1). This result is also 
useful for studying the multi-color solitons supported by this type of multistep parametric interaction (Kivshar, 
Alexander, and Saltiel [1999]; Huang [2001]). We note again the THG multistep cascading can be efficient only 
if the two processes can be simultaneously phase-matched. Due to dispersion in a bulk material, generally it is 
impossible to achieve double phase matching. However, due to the recent progress in the design of the QPM 
structures, the double phase matching can be achieved in nonlinear photonic materials (Zhang, Wei, Zhu, Wang, 
Zhu, and Ming [2001]; Luo, Zhu, He, Zhu, Wang, Liu, Zhang, and Ming [2001]; Fradkin-Kashi, Arie, Urcnski, and 
Rosenman [2002]), although there exist a number of technical problems to be solved. 



3.2. WAVELENGTH CONVERSION 

So far we discussed only the multistep cascading process that allows to generate the third harmonic by combining the 
SHG and SFG phase-matched parametric interactions in a single crystal. Let now consider another single-crystal 
multistep process that combines SHG and DFM and mimic in this way the four-wave mixing (FWM) process 
with two input waves lOp and Wsigi resulting in the generation of a signal at tJout ~ '2i^p ~ ^sig- The idea of this 
type multistep cascading is illustrated in Fig. 6. Usually, the difference Agig — Xp is smaller than 50 nm (with 
\p ~ 1550n7Ti) that leads to the result that when SHG process is phase matched the DFM process is very close 
to exact phase matching. Here we have the situation of double phase matched multistep cascading. This allows a 
very efficient conversion from Asig to Asig — Ap/2 which is in 10^ 10^ times larger than what can be obtained with 
the direct FWM process. To be more specific, let us consider the parametric interaction shown in Fig. 6 which is 
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described by the system of parametrically coupled equations, 



dz 

_ /l2„iAfcsHG2 A„ A A „iAfcDFGZ 

dz 



out 



dz 



(3.6) 



— - _ i(T5/i2/iQyje , 
dz 

where AfcsHG = A;2 ^ 2fci + G„i and A/sdfg = ^2 ^ ^"sig ^ ^'out , fi to (T5 are the couphng coefficients proportional to the 
second-order nonlinearity parameter dej / ■ The vector Gm is one of the QPM vectors used for achieving the phase 
matching (Fejer, Magel, Jundt, and Byer [1992]) in the QPM structure. If the birefringence phase matching is used, 
then Gm = 0. Two phase-matching parameters, i.e. A/cshg and A/cdfc, are involved into this parametric cascaded 
interaction. However, if the signal wavelength is sufficiently close to that of the pump, i.e. |Asig — Ap| <^ Ap, the 
tuning curves for the two processes practically overlap, and we have the situation in which if one of the parametric 
processes is phase matched, the other one is also phase matched. In other words, in this case the signal wavelength 
is sufficiently close to that of the pump, and we work under the conditions of double phase matching. If we neglect 
the depletion, the amplitude of the phase-matched SH wave can be found as ^2(2) = —i(J2ApZ. Then, the output 
signal can be written in the form, 

14 (T-llS ^ /l^2 2| /|2|2| . iz Sin"^ (AfcDFG-L/2) 

( AfcDFG j 



In the limit A/cdfg — * 0, the efficiency becomes 



lAoutiL^ ^ -ajal\Al\'\A,,,\^L\ (3.8) 



or 



Vout = \voPpL\ (3.9) 

where tjq is the normalized efficiency measured in W~^cm~^; it depends on the overlap integral between the 
interacting modes and the effective nonlinearity doff, and it is the same normalized value that describes the efficiency 
of the ffist step, the SHG process (Chou, Brener, Fejer, Chaban, and Christman [1999]). From Eq. (3.7) it follows 
that the output signal is a linear function of the input signal — the property important for communications. Also, 
the efficiency is proportional to the fourth power of the crystal length L. Detailed theoretical description of this 
cascading process can be found in Gallo, Assanto, and Stegenian [1997]; Gallo and Assanto [1999]; Chen, Xu, 
Zhou, and Tang [2002]. 

For the DFG process, the width of the phase-matching curve depends on the type of the crystal, its length, 
and the phase-matching method. Several proposal for increasing the phase-matching region have been suggested, 
including the use of the phase-reversed QPM structures (Chou, Parameswaran, Fejer, and Brener [1999]); the 
pump deviation from the exact phase matching (Chou, Brener, Parameswaran, and Fejer [1999]); the periodically 
chirped (phase modulated) QPM structures (Asobe, Tadanaga, Miyazawa, Nishida, and Suzuki [2003]); the phase 
shifting domain (Liu, Sun, and Kurz [2003]). In particular, the paper Gao, Yang, and Jin [2004] reports that 
the use of sinusoidally chirped QPM superlattices provides broader bandwidth and more flat response compared to 
homogeneous and segmented QPM structures. 

Novel cascaded x*-^'' wavelength conversion schemes are based on the SFG and DFG processes and the use of two 
pump beams, as proposed and demonstrated in Xu and Chen [2004]; Chen and Xu [2004]. The conversion efficiency 
is enhanced by 6 dB, as compared with the conventional cascaded SHG+DFG wavelength conversion configuration. 
The cascading steps in this SFG-I-DFG wavelength conversion method are u>sf = i^pi + Wp2, and Wout = ujsf — t^sig, 
so that the effective third-order interaction is the totally nondegenerate FWM process: Wout = ojpi + <jJp2 — ^sig. 

The second-order cascading wavelength conversion is one of the best examples of the richness of the multistep 
cascading phenomena. In the initial stage of development, this concept was experimentally demonstrated in the 
media employing two types of the phase matching techniques, i.e. the birefringence phase matching in a bulk crystal 
(e.g. Tan, Banfi, and Tomaselli [1993]; Banfi, Datta, Degiorgio, Donelli, Fortusini, and Sherwood [1998]) and 
periodically-poled nonlinear LiNbOa crystal (Banfi, Datta, Degiorgio, Donelli, Fortusini, and Sherwood [1998]). 
Since that, this parametric process was extensively investigated not only as an interesting multistep parametric 
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Figure 7: Right: periodically poled LiNbOa crystal for the cascaded wavelength conversion. Left: wavelength 
conversion with 1545 nm pump and four inputs in the range 1555-1560 nm (Chou, Brener, Lenz, Scotti, Chaban, 
Shmulovich, Philen, Kosinski, Parameswaran, and Fejer [2000]). 

effect that may occur in different nonlinear media (see Table 3), but also as a proposal for realistic all-optical 
communication devices (Chou. Brener, Fejer, Chaban, and Christman [1999]; Chou, Brener, Lenz, Scotti, Chaban, 
Shmulovich, Philen, Kosinski, Parameswaran, and Fejer [2000]; Kunimatsu, Xu, Pelusi, Wang, Kikuchi, Ito, and 
Suzuki [2000]). Nowadays, there is a real progress in the suggesting this device for optical communication industry 
(Cardakli, Gurkan, Havstad, Willner, Parameswaran, Fejer, and Brener [2002]; Cardakli, Sahin, Adamczyk, Willner, 
Parameswaran, and Fejer [2002]), because of clear advantages over all other devices used for the wavelength shifting. 
Indeed, one such device can simultaneously shift several channels. As is shown in Fig. 6, the spectra of the shifted 
signal is a mirror image of the origin. This feature can be used to invert the signal chirp for dispersion management 
in transmission systems. A successful experimental demonstration of this property has been reported recently in 
Kunimatsu, Xu, Pelusi, Wang, Kikuchi, Ito, and Suzuki [2000]: a 600-fs pulse transmission over 144 km using 
midway frequency inversion with this type second-order cascaded wavelength conversion resulted in a negligible 
pulse distortion. 

The waveguide made in the LiNbOa crystals (see Fig. 7) are by now the most suitable nonlinear structures for the 
cascaded simulation of the FWM wavelength shifting. The important features of this device is an almost perfect 
linear dependence between the input and output signals for more than 30 dB of the dynamic range, instanta- 
neous memory less transparent wavelength shifting that can be used at rates of several tera- Hertz, and transparent 
crosstalk- free operation (Cardakli, Sahin, Adamczyk, Willner, Parameswaran, and Fejer [2002]). Additionally, 
Couderc, Lago, Barthelemy, De Angelis, and Gringoli [2002] demonstrated that the wavelength conversion mul- 
tistep cascading system can support parametric solitons in the waveguiding regime, we discuss these results in 
more detail in Sec. 5.3 below. Experimental works reporting the second-order cascaded wavelength conversion are 
summarized in Table 3. 



Table 3: Experimental results on the multistep SHG and DFM cascading 



Nonlinear 
crystal 


Ap(As„„^) [mw] 


L [cm] 


Regime 


Phase 

matching 

method 


Refs. 


BBO 


1.064 (1.090) 


1 


Pulsed (30ps) 


BPM 


a 


MBA-NP 


1.064 (1.090) 


0.32 


Pulsed (30ps) 


BPM 


b 


LiNbOa 
waveguide 


1.533 (1.535) 


1 


CW&Pulsed 
(7ps) 


QPM 


c 


LiNbOs 


1.8 (1.863) 


1.9 


Pulsed (20 ps) 


QPM 


d 


NPP 


1.148 (1.158) 


0.28 


Pulsed (20 ps) 


NCPM 


e 


LiNbOa 
waveguide 


1.562 (1.600) 


4 


CW 


QPM 


f 


Ti:LiNb03 
waveguide 


1.103 (1.107) 


5.8 


Pulsed (20 ps) 


NCPM 


g, h 


LiNbOa 
waveguide 


1.545 (1.580) 


5 


CW 


QPM 


i 


Ti:LiNbOa 
waveguide 


1.556 (1.565) 


7.8 (8.6) 


CW & pulsed 
(6 ps) 


QPM 


j 
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Nonlinear 
crystal 




L [cm] 


Regime 


Phase 

matching 

method 


Refs. 


LiNbOa 
waveguide 


1.565 (1.585) 


2 


CW 


QPM 


k 


LiNbOg 
waveguide 


1.542 (1.562) 


1 


CW & pulsed 


QPM 


1 


LiNbOa 
waveguide 


1.550 (1.560) 




CW 


QPM 


m 


LiNbOa 
waveguide 


1.553 (1.565) 


6 


CW 


QPM 


n 


Ti:LiNb03 
waveguide 


1.557 (1.553) 


6 


CW 


QPM 


o 


LiNbOg 
waveguide 


1.532 (1.565) 


2 


CW 


QPM 


P 


LiNbOa 
waveguide 


1.537 (1570) 


4.5 


CW 


QPM 


q 


LiNbOa 
waveguide 


1.558.. 1.568 (1.600) 


3.4 


CW 


PM 
QPM^ 


r 


MgO:LiNb03 
waveguide 


1.543(1.573) 


5 


CW 


QPM 


s 


Ti:LiNbOa 
waveguide 


1.55(1.62) 


3 


CW 


QPM*^ 


t 


LiNbOa 
waveguide 


1.545(1.580) 


5 


CW 


QPM 


u 



References: 

" Tan, Banfi, and Tomaselli [1993] 

'' Nitti, Tan, Banfi, and Degiorgio [1994] 

Trevino-Palacios, Stegeman, Baldi, and De Micheli [1998] 

Banfi, Datta, Degiorgio, and Fortusini [1998] 
" Banfi, Datta, Degiorgio, Donelli, Fortusini, and Sherwood [1998] 
^ Chou, Brener, Fejer, Chaban, and Christman [1999] 
^ Cristiani, Banfi, Degiorgio, and Tartara [1999] 

Banfi, Christiani, and Degiorgio [2000] 
' Chou, Brener, Lenz, Scotti, Chaban, Shmulovich, Philen, Kosinski, Parameswaran, and Fejer [2000] 
-' Schreiber, Suche, Lee, Grundkotter, Quiring, Ricken, and Sohler [2001] 

* Cristiani, Liberale, Degiorgio, Tartarini, and Bassi [2001] 
' Ishizuki, Suhara, Fujimura, and Nishihara [2001] 

Cardakli, Sahin, Adamczyk, Willner, Parameswaran, and Fejer [2002] 
" Harel, Burkett, Lenz, Chaban, Parameswaran, Fejer, and Brener [2002] 
° Cristiani, Degiorgio, Socci, Carbone, and Romagnoli [2002] 
" Zeng, Chen, Chen, Xia, and Chen [2003] 
" Zhou, Xu, and Chen [2003] 

Asobe, Tadanaga, Miyazawa, Nishida, and Suzuki [2003] 
" Bracken and Xu [2003] 

* Gao, Yang, and Jin [2004] 
" Sun and Liu [2003] 



^ phase modulated QPM 

^three types of QPM grating are compared: homogeneous, segmented gratings and sinusoidally chirped 
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3.3. TWO-COLOR MULTISTEP CASCADING 

By two-color niultistcp cascading in a quadratic medium wc understand the multi phase-matched parametric inter- 
action between several waves which possess two frequencies (or wavelengths) only. One of the ways to introduce a 
parametric process involving more than one phase-matched interaction with two wavelengths is to consider a vec- 
torial interaction between the waves with different polarizations, or degenerate interaction between allowed modes 
in a waveguide. We denote two waves at the fundamental frequency (FF) (at A = Afund) as A and B, and the two 
waves of the second harmonic (SH) field (at Xsh = Xfund/'^ ), as S and T. Each pair of the eigenmodes [(A,B) and 
(S,T)] can be, for example, two orthogonal polarization states or two different waveguide modes at the fundamental 
and second-harmonic wavelengths, respectively. There exists a finite number of possible multistep parametric inter- 
actions that can coupled these waves. For example, if we consider the AA-S&AB-S cascading, then the multistep 
cascading is composed of the following sub-processes. First, the fundamental wave A generates the SH wave S 
via the type I SHG process. Then, by the down-conversion process SA-B, the other fundamental eigenmodc B 
is generated. At last, the initial FF wave A is reconstructed by the processes SB-A or AB-S, SA-A. When we 
deal with two orthogonal polarizations, the two principal second-order processes AA-S and AB-S are governed by 
two different components (or two different combinations of the components) of the x^^^ susceptibility tensor, thus 
introducing additional degrees of freedom into the parametric interaction. The classification of different types of 
the multistep parametric interactions has been introduced by Kivshar, Sukhorukov, and Saltiel [1999] and Saltiel, 
Koynov, Deyanova, and Kivshar [2000]. 



Table 4: Two-color multistep cascading processes 



Multistep- cascading 
schemes 


No of 
waves 


SHG processes 


Equivalent 

cascading 

schemes 


wc/sc 


Refs. 


AA-S: AB-S 


3 


Type I & 


I Type II 


BB-S:AB-S; 
AA-T:AB-T; 
BB-T:AB-T 


sc 


a,b,c 


AA-S:AB-T 


4 


Type I & 


c Type II 


BB-S:AB-T; 
AA-T:AB-S; 
BB-T:AB-S 


wc 


d,e 


AA-S:BB-S 


3 


Type I K 


s Type I 


AA-T:BB-T 


WC 


f,g,h,i 


AA-S:AA-T 


3 


Type I & 


s Type I 


BB-S:BB-T 


sc 


h,j 


AB-S:AB-T 


4 


Type II <( 


k Type II 




sc 




AA-S:BB-S:AB-S 


3 


Type I & 


J Type II 


AA-T:BB- 




k,l,m 



T:AB-T 



" Saltiel and Deyanova [1999] 

^ Saltiel, Koynov, Deyanova, and Kivshar [2000] 

" Petrov, Albert, Etchepare, and Saltiel [2001]; Petrov, Albert, Minkovski, Etchepare, and Saltiel [2002] 

DeRossi, Conti, and Assanto [1997] 
" Pasiskevicius, Holmgren, Wang, and Laurell [2002] 
•'' Assanto, Torelli, and Trillo [1994] 
" Kivshar, Sukhorukov, and Saltiel [1999] 

Grechin, Dmitriev, and Yur'ev [1999]; Grcchin and Dmitriev [2001b] 
* Grechin and Dmitriev [2001b] 

^ Trevino-Palacios, Stegeman, Demicheli, Baldi, Nouh, Ostrowsky, Delacourt, and Papuchon [1995] 

Trillo and Assanto [1994] 
' Towers, Sammut, Buryak, and Malomed [1999]; Towers, Buryak, Sammut, and Malomed [2000] 
Boardman and Xie [1997]; Boardman, Bontemps, and Xie [1998] 

Different types of the multistep parametric interactions can be divided into two major groups. The first group 
is composed by the parametric interactions with two common waves in both cascading processes, and the two 
processes are strongly coupled (SC). For the other group, both the parametric processes share one common wave, 
and these processes are weakly coupled (WC). The same classification can be applied to other types of the multistep 
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parametric interactions. In Table 4, we present an updated classification of two-color multistep parametric processes 
and the publications where they are analyzed. 

The first analysis of the multistep parametric interactions of this type has been carried out by Assanto, Torelli, 
and Trillo [1994] and Trillo and Assanto [1994]. These authors studied the simultaneous phase-matching of two 
SHG processes of the type AA-S and BB-S, where S denotes the SH wave whereas A and B stand for the the 
FF waves polarized in the perpendicular planes. The two orthogonal FF fields interact through the generated SH 
wave. All-optical operations and the polarization switching can be performed on the base of this scheme. Kivshar, 
Sukhorukov, and Saltiel [1999] demonstrated that this two-color parametric interaction can support two-color 
spatial solitary waves. Trevino-Palacios, Stegeman, Demicheli, Baldi, Nouh, Ostrowsky, Delacourt, and Papuchon 
[1995] considered the interference between the parametric processes AA-S and AA-T, where S and T are two different 
modes of the waveguide at the frequency 2a;. In both the cases mentioned above, two different parametric processes 
share one and the same fundamental wave. However, it is possible that the multistep cascading interaction involves 
the SHG process of the type AB-S. Such a four-wave multistep cascading process was considered by DeRossi, Conti, 
and Assanto [1997]. The two SHG processes were AB-S and AA-T, where (A,B) and (S,T) are pairs of the modes 
at the frequencies cu and 2u!, respectively. DeRossi, Conti, and Assanto [1997] concluded that a devise based on 
this type of multistep cascading can operate as an all-optical modulator or as an all-optical switch, with a good 
switching contrast at 1.55 /im. 

Another interesting process of the multistep interaction was considered by Boardman and Xie [1997]; Boardman, 
Bontemps, and Xie [1998] who studied the parametric mode coupling in a nonlinear waveguide placed in a magnetic 
field. In this case, the simultaneous coexistence of six SHG processes is possible, namely, ooo, ooe, oce, cec, ceo, and 
coo. The exact number of allowed parametric interactions depends on the symmetry point group of the material. 
It was shown that, by controlling a ratio of the input fundamental components, one of the SH components can be 
controlled and switched off. Importantly, repulsive and collapsing regimes for the interacting parametric solitons 
can be produced by switching the direction of the magnetic field. 

Saltiel and Deyanova [1999] considered the possibility to realize the efficient polarization switching in a quadratic 
crystal that supports simultaneous phase matching for both Type I and Type II SHG processes (e.g. ooo and oeo). 
At certain conditions, the fundamental beam involved in such a process can accumulate a large nonlinear phase 
shift at relatively low input power (Saltiel, Koynov, Deyanova, and Kivshar [2000]). The SHG process that can be 
realized by two possible pairs of the simultaneously phase-matched processes (ooo, ooe) and (ooe, eee) was studied 
theoretically in Grechin, Dmitriev, and Yur'ev [1999]. The effect of SHG by simultaneous phase-matching of three 
parametric processes (ooe, eee, and oee) in a crystal of LiNbOa was estimated theoretically in Grechin and Dmitriev 
[2001b]. It was shown that, at certain conditions, one can obtain polarization insensitive SHG. 

Also, we would like to mention several experimental studies of the two-color multistep cascading interactions. In 
the paper of Trevino-Palacios, Stegeman, Demicheli, Baldi, Nouh, Ostrowsky, Delacourt, and Papuchon [1995], an 
interplay of two Type I SHG processes with a common fundamental wave was observed. Petrov, Albert, Etchepare, 
and Saltiel [2001] and Petrov, Albert, Minkovski, Etchepare, and Saltiel [2002] performed the experiment with 
a BBO crystal in which, as a result of the simultaneous action of the Type I SHG and Type II SHG interactions, 
the generation of the wave orthogonally polarized to the input fundamental wave was observed. Couderc, Lago, 
Barthelemy, De Angelis, and Gringoli [2002] demonstrated that the multistep cascading interaction can support 
parametric solitons in the waveguiding regime. In this latter case, the multistep interaction simulates an effective 
nearly-degenerate FWM process and, in this sense, it is almost "two-color". Pasiskevicius, Holmgren, Wang, and 
Laurell [2002] realized experimentally the simultaneous SHG process that gives two SH waves with the orthogonal 
polarizations in the blue spectral region by use of the Type II and Type I QPM phase matching in a periodically 
poled KTP crystal. 

As a simple example of the two-color multistep parametric interaction, we consider the cascading of the Type 
I and Type II SHG processes according to the scheme: BB-S and AB-S. Here, the SH wave is generated by two 
interactions. Thus, we can expect an increase of the SHG efficiency when both A and B waves (i.e., the ordinary 
and extraordinary waves) arc involved. However, if only one of the waves, say the wave A, is launched at the 
input, this type of the double phase-matched interaction will lead to the generation of a wave perpendicular to the 
input wave, through the parametric process SA —f B. Additionally, the fundamental wave A accumulates a strong 
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nonlinear phase shift. This parametric interaction is described by the foUowing equations for plane waves, 

dz 

^ = -in^A^e'^^^^^'- - iaiABe'^^^^^\ (3.10) 
dz 

dz 

where A, S, and B are the complex amplitudes of the input fundamental wave, the second-harmonic wave, and the 
orthogonally polarized wave at the fundamental frequency, cri and (T2 are defined above, and 

27r(icft,ii / ^^1,2,1 

f3,4,5 ~ 



If we neglect dispersion of the index of refraction, i.e. ua — ns — ^^2, then we can assume that ai ~ (T2 and ~ 
CTs ~ cr4/2. The phase-mismatch parameters are defined as AfcgHG ~ fc2 + and AfcopQ ~ k2 — kA — kB+Gq, 
where Gp and Gq are two QPM vectors used for the phase matching. 

Solution of this system, with respect of the amplitude of the wave B, has been obtained by Saltiel and Deyanova 
[1999]) in the approximation of nondepleted pump. It gives the following result, 

= nZ' \A\'Asm{D^L)e^^-\ (3.11) 

where ^ 

D± = ^(AfcsHG - AfcDFo) ± -XT^l^f- 

2 AfcsHG 

Equation (3.11) shows that the generation of the component B by this multistcp interaction mimic the FWM 
process of the type AAA* — governed by the cascaded cubic nonlinearity. This effective cubic nonlinearity 
leads to the accumulation of a nonlinear phase shift by the wave A that, similar to the case of the cascaded THG 
process, includes a contribution of high-order (> 3) nonlinearities (Saltiel, Koynov, Deyanova, and Kivshar [2000]). 
Existence of two-color multistcp parametric solitons and waveguiding effects are the result of the nonlinear phase 
shift collected by the interacting waves (Kivshar, Sukhorukov, and Saltiel [1999]). 



3.4. FOURTH-HARMONIC MULTISTEP CASCADING 

Another fascinating example of the multistcp cascading interaction is the fourth-harmonic generation (FHG) in 
a single crystal with the second-order nonlinearity. There are two possible parametric processes of the second 
order that lead to the generation of a fourth harmonic. In both the cases the fourth harmonic (FH) amplitude is 
proportional to the factor (x^^'')'^: 

1. SHG + SHG: lu + lo = 2uj: 2lu + 2lu = Auj ; 

2. SHG + SFG -l-SFG: w -f w = 2tj; 2tj -f w = 3w; 3w -t- w = 4w. 

The phase matching conditions define which of these processes will be more effective. Obviously, the former process 
is easier to realize technically because it requires only two phase-matching conditions to be fulfilled simultaneously, 
and therefore many papers deal with this case. 

The first experiment on the generation of a forth-harmonic wave by cascading was reported in the paper by 
Akhmanov, Dubovik, Saltiel, Tomov, and Tunkin [1974], where the FHG cascaded process in Lithium Formiate 
crystal was studied. As pointed out by the authors, the generation is due to the simultaneous action of two and 
three processes with the involvement of the quadratic and cubic nonlinearities. The FHG process in CdGeAs2 
was observed by Kildal and Iselcr [1979]. Both these pioneer papers were motivated by the idea to estimate the 
magnitude of the direct fourth-order nonlinearity in terms of (x*-^^)'^ . One of the first studies that reported the 
efficiency of the cascaded FHG process is the paper by Hooper, Gauthier, and Madey [1994] where the efficiency of 
3.3 X 10~^% was obtained in a single LiNbOa crystal. The other paper by Sundheimer, Villencuve, Stegeman, and 
Bierlein [1994b] reported the efficiency of 0.012%. However, in those experiments only the first step uj + uj = 2uj 
was phase-matched while the other one was not matched, leading to the overall low conversion efficiency. Somewhat 
larger efficiency of 0.066% was reported in the work of Baldi, Trcvino-Palacios, Stegeman, Demicheli, Ostrowsky, 
Delacomt, and Papuchon [1995] where a periodically poled LiNbOa waveguide was used. Both steps of the multistep 
parametric interaction, namely uj + uj ~ 2uj and 2uj + 2lu = 4a;, were phase matched; the first step was realized 
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through the first-order QPM process, while the second step was realized through the 7-th order QPM process. 
As is shown by Norton and de Sterke [2003b] and Sukhorukov, Alexander, Kivshar, and Saltiel [2001], if both 
steps are phase-matched, the resulting efficiency should be close to 100%. In the second paper, the existence and 
stability of the normal modes for such a multistep cascading system have been studied. Some possibilities for the 
double-phase-matched FHG process for certain input wavelengths in single-crystals of LiNbOs, LiTaOa, KTP, and 
GaAs have been shown by Pfistcr, Wells, Hollberg, Zink, Van Baak, Lcvenson, and Bosenbcrg [1997] and Grcchin 
and Dmitriev [2001a]. The possibility for FHG by double phase-matching in broader spectral region by use of the 
phase-reversed QPM structure was discussed by Sukhorukov, Alexander, Kivshar, and Saltiel [2001]. 

In a very interesting experiment, Broderick, Bratfalean, Monro, Richardson, and dc Sterke [2002] demonstrated 
cascaded FHG in a 2D nonlinear photonic crystal with the efficiency 0.01% in a (not optimized) 2D planar QPM 
structure. Several useful efficient schemes for FHG in 2D nonlinear photonic crystals have been proposed by Saltiel 
and Kivshar [2000a] and dc Sterke, Saltiel, and Kivshar [2001]. An optimal design of 2D nonlinear photonic 
crystals for achieving the maximum efficiency of FHG has been discussed by Norton and dc Sterke [2003a] and 
Norton and dc Sterke [2003b]. The use of the FHG multistep interaction for the frequency division schemes (4:1) 
and (4:2) have been discussed in Dmitriev and Grechin [1998] and Sukhorukov, Alexander, Kivshar, and Saltiel 
[2001]. 

The basic equations describing the FHG parametric process in a double-phase-matched QPM structure for the 
interaction of plane waves can be written in the form (see, e.g.. Hooper, Gauthicr, and Madey [1994]) 

dz 

. .2 iAfcsHGZ ■ A A* -iAkiZ (o l ON 

dz 

—j— — —1(77/126 , 

where Ai, A2 and A4 are the complex amplitudes of the fundamental, second- harmonic, and fourth-harmonic waves, 
respectively. The parameters ci and a2 are defined above, and 

47r 

CTe.y — -: «off,ii. 

Mn2A 

As before, if we neglect the index of refraction dispersion (ni ~ n,2 ~ 71,4) then we can accept that ci ~ <J2 and 
(Tg ~ (T7. Phase mismatch parameters are AfcsHG = ^2 — 2fcA -I- Gp and Ak^ = k^ — 2k2 + Gq , where Gp and Gq are 
two QPM reciprocal vectors. The role of high order nonlinearities are not included in (3.12) since their contribution 
is rather small when one works in conditions close to double or triple phase-matching. 

Solution of the system (3.12) can be found neglecting the depletion effects. It reveals that the phase-matched 
FHG wave is generated when either one of the following phase-matching conditions is satisfied, 

1. AfcsHG 0; 

2. Afc4 — > 0; 

3. AfcsHG + ^ki-k2- 2ki +Gp + Gq — > 0; and 

4. 2AfcsHG + Afc4 = k4- 4/ci + 2Gp + Gq — >0 

The case AfcgHc — > gives the strongest phase-matching process among all those cases, and the generated FH wave 
exceeds by several orders of magnitude FH generated by other schemes. For the double-phase-matching process 
(when AfcgHG — * and Afc4 — > 0) the squared amplitude of forth harmonic is given by the expression (de Sterke, 
Saltiel, and Kivshar [2001]): 

\A,{L)\^^laWMi\'L'. (3.13) 

Introducing again the normalized efficiency (measured in W^^cm"^) for the ffist and second steps as 770,1 and 
ilo,2, respectively, we can obtain the efficiency of the cascaded FHG process, 

= lvlim,2PlL''- (3.14) 

Thus, the efficiency of the cascaded FHG process in a single x^^"* is proportional to the 6-th power of the length 
and the 3-rd power of the pump, and it can be estimated with the known efficiencies for the separated steps. 
For the pump intensities at which the depletion effect of the fundamental and second-harmonic waves can not be 
neglected, the system (3.12) should be solved numerically. As shown in Zhang, Zhu, Zliu, and Ming [2001], the 
100% conversion of the fundamental wave into the fourth- harmonic wave is possible independently on the ratio of 
the nonlinear coupling coefficients 02 1 o^. This behavior is in contrast with the x^^^ -based cascaded THG process 
where the 100% conversion is possible only for a specific ratio of the nonlinear coupling coefficients (see Sec. 3.1). 
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Figure 8: The timing cliaracteristics of OPO with the simuhancously phase matched SHG and SFG processes, (a) 
SFG between the signal and pump, and (b) SHG of the signal. Inset: output power vs. wavelength (Zhang, Hebling, 
Kuhl, Ruble, Palfalvi, and Giessen [2002]). 



3.5. OPO AND OPA MULTISTEP PARAMETRIC PROCESSES 

Many studies dealing with optical parametric oscillators (OPOs) and optical parametric amplifiers (OPAs) have 
observed, in addition to the expected signal and idler waves, other waves at different wavelengths that are coherent 
and are generated with a good efficiency. These waves were identified as being a result of additional phase-matched 
(or nearly phase-matched) second-order parametric processes in OPOs and OPAs. In some of the cases, up to five 
additional waves coming out from the nonlinear crystal were detected and measured. By tuning the wavelength of 
the signal and idler waves simultaneously, the additional outputs from OPO can be frequency tuned as well. An 
example is shown in Fig. 8 adapted from Zhang, Hebling, Kuhl, Ruble, Palfalvi, and Giessen [2002]. The study of 
the simultaneous phase matching in OPO was started more than 30 years ago. As shown in Ammann, Yarborough, 
and Falk [1971], in many nonlinear crystals the phase-matching tuning curves for OPO and SHG of the signal and 
the idler cross and the experimental observations of the corresponding double phase-matching parametric processes 
are possible. 

Plane-wave single-pass theory of the new frequency generation with OPO with the simultaneous frequency dou- 
bling of the signal wave in the OPO crystal was presented in Aytur and Dikmelik [1998]. The case of OPO with the 
simultaneous frequency mixing between the pump and signal waves — in Dikmelik, Akgun, and Aytur [1999] and 
Morozov and Chirkin [2003]. Moore, Koch, Dearborn, and Vaidyanathan [1998] considered theoretically the simul- 
taneous phased matched tandem of OPOs in a single nonlinear crystal. The signal wave of the first OPO process 
becomes a pump wave of the second process. The simultaneous action of the parametric generation, uop uoi + ujg, 
and SFG between the pump and signal wave, cjsf = + Ws, was also analyzed by Huang, Zhu, Zhu, and Ming 
[2002] , where the main goal was to show that this type of the multistep cascading can be used for the simultaneous 
generation of three fundamental colors. 

Table 5: Experimental results on the OPO and OPA multistep cascading 

Nonlinear Phase matching Phase L Regime Refs 
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method 
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Nonlinear 
crystal 



Simultaneously 
PM processes 



Phase 

matching 
method 



Regime Refs 



BBO 
KTP 

LiNbOs 



LiNbOS 



/3-BBO 



KTP 



KTP 



LiNbOs 



KTA 



0.74-0.89 
0.82-0.92 

1.064 



LiNbOs 0.532 
LiNbOs 0.78-0.8 



0.793 



0.53 



0.74-0.76 



LiNbOs 1.064 
LiNbOs 0.79-0.81 
KTP 0.827 
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LiTaOs 0.532 
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1.064 
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" Andrews, Rabin, and Tang [1970] 
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Ammann, Yarborough, and Falk [1971] 

Bakker, Planken, Kuipers, and Lagendijk [1989] 

Powers, Ellingson, Pelouch, and Tang [1993] 
" Schiller and Byer [1993] 
■'' Petrov and Noack [1995] 

Hebling, Mayer, Kuhl, and Szipocs [1995] 
^ Myers, Eckardt, Fejer, Byer, Bosenberg, and Pierce [1995] 
' Schiller, Breitenbach, Paschotta, and Mlynek [1996] 

Butterworth, Smith, and Hanna [1997] 

Burr, Tang, Arbore, and Fejer [1997] 
' Varanavicius, Dubietis, Berzanskis, Danielius, and Piskarskas [1997] 

Kartaloglu, Koprulu, and Aytur [1997] 
" Vaidyanathan, Eckardt, Dominic, Myers, and Grayson [1997] 
° McGowan, Reid, Penman, Ebrahimzadeh, Sibbett, and Jundt [1998] 
^ Koprulu, Kartaloglu, Dikmelik, and Aytur [1999] 
« Zhang, Hebling, Kuhl, Ruble, and Giessen [2001] 

Zhang, Hebhng, Kuhl, Ruble, Palfalvi, and Giessen [2002] 
" Du, Zhu, Zhu, Xu, Zhang, Chen, Liu, Ming, Zhang, Zhang, and Zhang [2002] 
' Zhang, Hebling, Bartels, Nau, Kuhl, Ruble, and Giessen [2002] 
" Xu, Liang, Li, Yao, Lin, Cui, and Wu [2002] 
" Kartaloglu and Aytur [2003] 
™ Kartaloglu, Figen, and Aytur [2003] 

Lee, Zhang, Huang, and Pan [2003] 

Abbreviations: BPM - birefringence phase matching; NCPM - noncritical phase matched; QPM - uniformly poled 
quasi-phase-matched structure; QPOS - quasi-periodical optical superlattices; APQPM - aperiodically poled QPM 
structure. 



In Table 5, we have summarized different experimental results on the multistep cascading processes observed in 
OPOs and OPAs. For each of this work, we also mention the additional second-order parametric processes. 

In the special case when all frequencies become phase related, i.e. : (jJ^ : cj^ = 3 : 2 : 1, OPO displays the 
unique properties (Kobayashi and Torizuka [2000]; Longhi [2001b]). As a matter of fact, this case corresponds 
to the third-harmonic multistep cascading (see Sec. 3.1) but realized in a cavity. The second harmonic of the idler 
simultaneously generated in the OPO crystal (or by externally frequency doubling) L02i ~ uji+uji will interfere with 
the signal wave of the frequency LUi. The resulting beat signal can be used for locking OPO at this particular tuning 
point and for realizing the frequency divisions (3:2) and (3:1), or vice versa. More details about the frequency 
division with these types OPO can be found in the papers by Lee, Klein, Meyn, Wallenstein, Gross, and Boiler 
[2003]; Douillet, Zondy, Santarelli, Makdissi, and Clairon [2001]; Slyusarev, Ikegami, and Ohshima [1999]; Zondy, 
Douiliet, Tallet, Ressayre, and Le Berrc [2001] ; Zondy [2003] . For the double phase-matched parametric interactions 
for the frequencies 3a;; 2a;; w, the total conversion from a 3a; wave to a 2a; wave was predicted when the pump is 
the 3a; wave (Komissarova and Sukhorukov [1993]), and from a 2a; wave to a 3a; wave, when the pump is the 2a; 
wave (Volkov and Chirkin [1998]) for the OPA scheme. Additionally, the theoretical studies in Longhi [2001b, a, c] 
show that OPO, when the signal and idler wave have the frequencies uj and 2a;, respectively, can produce different 
types of interesting patterns including the spiral and hexagonal patterns. The internally pumped OPOs, due to the 
simultaneous action of the SHG process and the parametric down-conversion, also shows the formation of spatial 
patterns and the parametric instability dynamics (Lodahl and Saffman [1999]; Lodahl, Bache, and Saffman [2000, 
2001]). 

Additionally, the multistep cascaded OPOs and OPAs can be an efficient tool for generating the entangled and 
squeezed photon states. In particular, Smithers and Lu [1974] considered the quantum properties of light generated 
during simultaneous action of the parametric processes LUp LOi + ujs and a;i = un+tOp. The theoretical predictions 
by Marte [1995b, a]; Eschmann and Marte [1997] suggest that, due to the multistep cascading in internally pumped 
OPO, this type of OPO can be an excellent source for generating twin photons with the sub-Poissonian statistics, 
and the generated SH wave exhibits a perfect noise reduction. Several other phase-matched schemes for generating 
the entangled and squeezed photon states have been proposed by Chirkin [2002]; Nikandrov and Chirkin [2002b] 
and Chirkin and Nikandrov [2003] who utilizes OPA with two simultaneously phase- matched processes, 2a; = a; + a; 
and 3a; = 2a; -f-a;. In the theoretical study (Nikandrov and Chirkin [2002a]), a possibility of generating squeezed 
light with the utilization of a single quadratic crystal was compared with that of two different crystals for each step. 
The conclusion is that a single-crystal cascaded method is more efficient. New possibilities for the generation of 
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Figure 9: Left: the phase matching curves for three simultaneously phase-matched processes in a single aperiodically- 
poled LiTaOa crystal. Right: visible red, green and blue outputs from a nonlinear crystal diffracted by a prism 
(Liao, He, Liu, Wang, Zhu, Zhu, and Ming [2003]). 

squeezed polarized light have been discussed by Dmitriev and Singh [2003] who considered the five-wave generation 
via the four simultaneously phase-matched parametric processes. 

3.6. OTHER TYPES OF MULTISTEP INTERACTIONS 

The multistep parametric interactions allow building compact frequency converters with several visible beams as 
the output. He, Liao, Liu, Du, Xu, Wang, Zhu, Zhu, and Ming [2003] reports on the simultaneous generation of 
all three " traffic signal lights" . The simultaneous generation of a pair blue and green waves have been achieved in 
Capmany, Bermudez, Callejo, Sole, and Dieguez [2000] by exploring self-doubling and self-frequency mixing active 
media. 



Table 6: Experiments on the simultaneous generation of several visible 
beams 
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\pm] 


COLORS 


PM 

method 


L [cm] 


Regime 


Ref. 


KTP 


1.023 


RED, 




0.45 


CW 


a 


waveguide 


0.716 


GREEN, 


BPM/QPM 












BLUE 










KTP 


1.650 


RED, 


QPM 


0.35 


pulsed 


b 


waveguide 




GREEN, 






(6 ps) 








BLUE 










LiNbOg 


1.620 


RED, 


QPM 




pulsed 


c 


waveguide 




GREEN, 






(7 ps) 








BLUE 










NYAB 


1.338 


RED, 


BPM 


0.5 


CW 


d 




0.807 


GREEN, 












0.755 


BLUE 












1.084 


GREEN, 


APQPM 


0.095 


CW 


e 


NdiLiNbOg 


0.744 


BLUE 












1.372 


RED, OR- 


APQPM 


0.3 


CW 


f 


NdiLiNbOs 


1.084 


ANGE, 












0.744 


GREEN, 














BLUE(2) 










SBN 


1.34 


RED, 


APQPM 


0.7 


CW 


g 




0.88 


GREEN, 














BLUE 










LiTaOg 


1.342 


RED, 


APQPM 


1 


CW 


h 




1.064 


GREEN, 














BLUE 
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Nonlinear \p COLORS PM L [cm] Regime Ref. 

crystal [a"^] method 

LiTaOg 1.342 RED, YEL- APQPM 1 CW i 

1.064 LOW, 
GREEN 



" Laurell, Brown, and Bierlein [1993] 
Sundheimer, Villeneuve, Stegeman, and Bierlein [1994b] 

Baldi, Trevino-Palacios, Stegeman, Demicheli, Ostrowsky, Delacourt, and Papuchon [1995] 

Jaque, Capmany, and Garcia Sole [1999] 
" Capmany, Bermudez, Callejo, Sole, and Dieguez [2000] 
•'' Capmany [2001] 

^ Romero, Jaque, Sole, and Kaminskii [2002] 

Liao, He, Liu, Wang, Zhu, Zhu, and Ming [2003] 
' He, Liao, Liu, Du, Xu, Wang, Zhu, Zhu, and Ming [2003] 

Recently, several studies presented successful attempts to obtain the simultaneous generation of red, green and 
blue radiation(the so-called RGB radiation) from a single nonlinear quadratic crystal. This is an important target for 
building compact laser-based projection displays. Theoretically, the parametric process for achieving the generation 
of three primary colors as OPO outputs was considered in Huang, Zhu, Zhu, and Ming [2002]. Liao, He, Liu, Wang, 
Zhu, Zhu, and Ming [2003] used a single aperiodically-poled LiTaOa crystal for generating 671, 532 and 447 nm 
(see Fig. 9) with three simultaneously phase-matched processes: SHG of 1342 and 1064 nm (the output of the dual- 
output Nd:YV04 laser) and SFG of 671 and 1342 nm. Jaque, Capmany, and Garcia Sole [1999] realized a different 
method for achieving the three-wavelength output. In their experiments, a nonlinear medium is the NdiYAla (603)4 
crystal that was pumped by a Ti:sapphire laser. The red signal at 669 nm was obtained by self-frequency doubling of 
the fundamental laser line. The green signal at 505 nm and a blue signal at 481 nm were obtained by self-SFG of the 
fundamental laser radiation at 1338 nm and the pump radiation (807 nm, for green, and 755 nm, for blue). All three 
processes were simultaneously phase matched by birefringence phase matching due to an exceptional situation that 
the three phase-matchings appear extremely close to each other and their tuning curves overlap. The generation 
of red, green and blue signals by triple phase-matching in LiNbOa and KTP periodically poled waveguides was 
reported also by Sundheimer, Villeneuve, Stegeman, and Bierlein [1994a] and Baldi, Trevino-Palacios, Stegeman, 
Demicheli, Ostrowsky, Delacourt, and Papuchon [1995]. The experimental efforts to build optical devices with the 
simultaneous generation of several visible harmonics arc summarized in Table 6. 

3.7. MEASUREMENT OF THE x^^^-TENSOR COMPONENTS 

As an important application of the multistcp parametric processes in nonlinear optics, we would like to mention the 
possibility making a calibration link between the second- and third-order nonlincarity in a nonlinear medium. Both 
THG and FWM processes in non-ccntrosymmetric nonlinear media can be used for this purpose. Measurements 
can be done in the phase-matched or non-phase-matched regimes. The non-phase-matched regime allows achieving 
a higher accuracy, however, it is more complicated since the signal is weak and additional care should to be taken 
to avoid the influence of the respective third-order effect in air. The basic idea of these types of measurements is 
to compare the THG or FWM signal in several configurations that include a proper choice of the direction and 
polarization of the input and output waves. In some of the configurations the output signal is generated in result 
of the direct contribution of the inherent cubic nonlincarity of a sample. In other cases, the signal is generated 
due to the cascade contribution, while in the third group, the contribution of both direct and cascaded processes is 
comparable. In this way, we can access the ratio 

x(2)(-3w,2w,w)x(2)(-2w,w,w) 
for the case of the THG multistcp interaction, and the ratio 
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for the case of the FWM multistep cascading. Because the information about the x*-^-* components is more available, 
we can determine quite accurate parameters of the cubic nonlinearity of non-centrosymmetric materials by using 
this internal calibration procedure which does not require the knowledge of the parameters of the laser beam. 

To illustrate that let us consider the phase-matched THG. In condition of non-depletion of the fundamental and 
second harmonic wave and neglecting the temporal and spatial walk off effect the equations (3.1) has following 
solution for phase matching A/cthg = ^3 ~ 3/ci — > 0: 



Apparent cubic nonlinearity consists two parts direct and cascading: 



(3) 
Xtot 



Aeff.dir 



(3) 
'^cff .case ' 



(3.15) 



(3.16) 



with 



(3) 
'^eff.casc 



WTTdcsjd, 



■cffjl 



AiniAfcsFG 



(3) 

One of the ways to separate the contribution of the two nonlincarities and express of x^g dh- terms of the 
product ((ieff,i)('^cff,ii) is to use the azimuthal (Banks, Feit, and Perry [2002]) or input polarization (Kim and 



Yoon [2002]) dependence of xlff dir ^^'^ case- "^^^ other way is to compare the TH signal obtained under the 



.(3) 



condition A/cthg 



with the TH signal obtained under one of the conditions AfcsHC — > or A/csfg 



0, 



where the signal is only proportional to the factor [xla casd^ (Akhmanov. Meisncr. Parinov, Saltiel, and Timkin 



(3) 

[1977]; Chemla, Begley, and Byer [1974]), and calculate Xeff dir' however this procedure gives two possible values 
due to an indeterminate sign. Considering all symmetry classes Fevc, Boulanger, and Douady [2002] found the 
crystal directions for which the second-order cascade processes give no contribution and, therefore, they are suitable 

(3) 

for the measurement of the value of Xcff dir- 

The parametric interaction that occurs when the degenerated FWM (DFWM) process is in non-centrosymmetric 
media is special because it consists of the steps of the optical rectification and linear electro-optic effects (Bosshard, 
Spreiter, Zgonik, and Gunter [1995]; Unsbo [1995]; Zgonik and Guntcr [1996]; Biaggio [1999, 2001]). Then, 
the measured x'^' component can be expressed through the squared electro-optic coefficient of the medium. The 
cascaded x'''^^ contribution in the crystals with a large electro-optic effect leads to a very strong cascaded DFWM 
effect which can exceed by many times the contribution of the inherent direct x*"^"* nonlinearity (Bosshard, Biaggio, 
St, FoUonier, and Gunter [1999]). 

Table 7 presents a summary of the experimental results on the measurement of the cubic nonlincarities by the 
use of the cascaded THG or cascaded FWM processes. 



Table 7: Experimental results for the x'^'-tensor components measured 
through the second-order multistep cascading processes 



Nonlinear 


^fund 


PM/NPM/ 


Cascading 


Reference 


crystal 


[fim] 


NCPM* 


scheme 




ADP 


1.06 


PM 


THG 


Wang and Baardsen [1969] 


GaAs 


10.6 


NPM 


FWM 


Yablonovitch, Flytzanis, and 
Bloembergen [1972] 


CdGeAs2 


10.6 


PM 


THG 


Chemla, Begley, and Byer [1974] 


KDP 


1.064 


PM 


THG 


Akhmanov, Meisner, Parinov, 
Saltiel, and Tunkin [1977] 


Qf-quartz 


1.91 


NPM 


THG 


Meredith [1981] 


/3-BBO 


1.054 


PM 


THG 


Qiu and Penzkofer [1988] 


/3-BBO 


1.053 


PM 


THG 


Tomov, Van Wonterghem, and 
Rentzepis [1992] 


KTP 


1.06 


PM 


DFWM 


DeSalvo, Hagan, Sheik-Bahae, 
Stegeman, Vanstryland, and 
Vanherzeele [1992] 


/3-BBO KD*P 


1.055 


PM 


THG 


Banks, Feit, and Perry [1999, 2002] 



d-LAP 
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Nonlinear \fund PM/NPM/ Cascading Reference 

crystal [m^^] NCPM* scheme 



KTP 



DAST# 

a-quartz 

KNbOa 

KTaOa 

SF59 BK7 
fused silica 

KNbOs 
DAST 
BaTiOs 
AANP* 



1.62 



1.064 

1.064 
1.318 
1.907 
2.100 

1.06 



1.390 
1.402 



NCPM 

PM 
NPM 



KDP 
LiNbOs 



BBO 1.064 
0.532 



PM 



NCPM 



NPM 



THG Feve, Boulanger, and Guillien 

[2000]; Boulanger, Feve, Delarue, 
Rousseau, and Marnier [1999] 

DFWM Bosshard, Biaggio, St, FoUonier, 
and Gunter [1999] 

THG Bosshard, Gubler, Kaatz, Mazer- 

ant, and Meier [2000] 



DFWM Biaggio [1999, 2001] 



FWM 



DFWM 



Taima, Koinatsu, Kaino, 
Franceschina, Tartara, Banfi, 
and Degiorgio [2003] 

Ganoov, Kulagin, Ryasnyanskii, 
Tugushev, and Usmanov [2003] 



* PM - phase matched; NPM - non phase-matched; NCPM - noncritically phase-matched; 

* organic crystals. 

§ 4. Phase-matching for multistep cascading 

In general, the simultaneous phase-matching of several parametric processes is hard to achieve by using the tradi- 
tional phase-matching methods, such as those based on the optical birefringence effect, except some special cases 
discussed above. However, the situation becomes quite different for the nonlinear media with a periodic variation of 
the sign of the quadratic nonlinearity, as it occurs in the fabricated one-dimensional (ID) QPM structures (Fejer, 
Magel, Jundt, and Byer [1992]) or two-dimensional (2D) x*-^^ nonlinear photonic crystals (Berger [1998]; Broder- 
ick, Ross, Offerhaus, Richardson, and Hanna [2000]; Salticl and Kivshar [2000a]). In this part of the review, we 
describe the basic principles of the simultaneous phase-matching of two (or more) parametric processes in different 
types of ID and 2D nonlinear optical supcrlattices. 

4.1. UNIFORM QPM STRUCTURES 

A bulk homogeneous nonlinear crystal possessing a quadratic nonlinearity is usually homogeneous everywhere. 
Several methods have been suggested and employed (Fejer [1998]) for creating a periodic change of the sign of the 
second-order nonlinear susceptibility 0?*^^^ as occurs in the QPM structure shown in Fig. 10. From the mathematical 
point of view, such a periodic sequence of two domains can be described by a simple periodic function, 

d{z) ^doY, 5rne'^'"^ (4.1) 

g,n = f — ) sm{TTmD), (4.2) 
\rmr / 

where Gm ~ (27rm/A) is the reciprocal QPM vector. The uniform QPM structure is characterized by a set of the 
reciprocal vectors: ±27r/A, ±47r/A, ±67r/A, ±87r/A, which can be used to achieve the phase- matching conditions 
when Ak — > 0. The integer number m (that can be both positive and negative) is called the order of the wave 
vector phase-matching. According to Eq. (4.1), the smaller is the order of the QPM reciprocal wave vector, the 
larger is the effective nonlinearity. If the filling factor D = 0.5, the effective quadratic nonlinearities (proportional 
to the parameter do5m) that correspond to the even orders QPM vectors vanish. Importantly, such uniform QPM 
structures can be used for simultaneous phase-matching of two parametric processes when the interacting waves 
are collinear or non-coUincar to the reciprocal wave vectors of the QPM structure. 
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h h 



D= li/A 



7 



d^'^ < 



Figure 10: Schematic of the uniform QPM structure: G„ 
factor, and A is the period of the QPM structure. 



is the reciprocal vector of the m-th order; D is the filhng 



4.1T. Colhnear case with two commensurable periods 

As an example of the multistep QPM interaction, we consider the third-harmonic multistep-cascading process 
under the condition that the interacting waves are collinear to the reciprocal wave vectors of the QPM structure. 
We denote the mismatches of the nonlinear material without modulation of the quadratic nonlincarity ("bulk 
mismatches") as A&i and A&2, where A6i — k2 ~ 2ki and A&2 = fcs — ^2 ~ ^i, and choose the period of the 



QPM structure in order to satisfy the phase-matching conditions Gmi — — A&i and = — A&2 

parametric processes, SHG and SFG, are characterized by the wave-vector mismatches Afci = k2 — 2fci 



The two 
Gmi and 



Afc2 = fca — fc2 — fci + Gm2! respectively, and they become simultaneously phase-matched for this particular choice 
of the QPM period. A drawback of this method is that it can satisfy simultaneously two phase-matching conditions 
for discrete values of the optical wavelength only. In particular, the values of the fundamental wavelength A for the 
double phase-matching condition can be found from the relation 



A52/TO2 - A6i/mi = 0, 



(4.3) 



where both Abi and A62 are functions of the wavelength. Equation (4.3) is valid for any pair of the second-order 
parametric processes. For the third-harmonic multistep-cascading process, Eq. (4.3) is transformed to the following, 



mi [3n(3a;) — 27i(2a;) — n{Lu)] — Irrii [7'i(2cij) — ?i(ti')] = 0, 



(4.4) 



where the arguments shows the wavelength dependence of the refractive index. For a chosen pair of integer numbers 
(mi,TO2), the required QPM period is found from the relation A = 27r|mi/A6i| or A = 27r|m2/A&2|. Such a method 
was used for double phase-matching by many researchers for the cascaded single-crystal third- and fourth-harmonic 
generation (see details and references in Table 2). The highest efficiency achieved with the uniform QPM structure 
so far is 19.3%, and it is achieved by the use of the combination of SHG (lst-order)/SFG(3rd-ordcr) and is reported 
in He, Liu, Luo, Jia, Du, Guo, and Zhu [2002]. 



4.1.2. Non-collincar case 

Collinearity between the optical waves and the reciprocal vectors of the QPM structure is an important requirement 
for achieving a good overlapping of all the beams and a good conversion efficiency. However, phase-matching is 
possible even in the case when some of the waves propagate under a certain angle to the direction of the reciprocal 
vectors of the QPM structure, i.e. for the non-coUinear case. With this method, the double phase-matching 
processes can be realized in a broad spectral range (Saltiel and Kivshar [2000b]). Such a type of non-collincar 
interaction will be efficient for the distances corresponding to the overlap of the interacting beams. 

As an example, we consider the THG multistep parametric process in the noncollinear geometry as shown in 
Fig. 11. in this case, the simultaneous phase-matching of two parametric processes lo -\- = 2io and a; -|- 2a; = 3w 
is required. We assume the fundamental wave at the input. As shown in Fig. 11, for the first process the phase- 
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matching can be achieved by the use of the the reciprocal vector G,„j , and the generated SH wave with the 
wavevector k2 is not colhnear to the fundamental wave: 2ki — G„jj = k2. The second process can be phase- 
matched by using the vector Gmj! ki + k2 — Gm2 = ^3. From Fig. 11, we can derive a result for the period of 
the QPM structure that allows achieving the double phase-matching of the processes of cascaded THG in a single 
QPM structure, 



2mi(fcg - 9kl) - 3(mi + m2)(fci - Akf) 

A = 27r ■ — . (4.5j 

mi(mi + m2)[2m2 — mi) 

The conditions for the double phase-matching for the THG multistep-cascading process in LiTaOa for the case when 
all waves are polarized along the z axis of the crystal were considered by Saltiel and Kivshar [2000b] . 



1/2 



4.2. NON-UNIFORM QPM STRUCTURES 

Non-uniform QPM structures also allow the simultaneous phase-matching of two parametric nonlinear processes. We 
consider three types of such non-uniform QPM structures: (i) phase-reversed QPM structures (Chou, Parameswaran, 
Fejer, and Brener [1999]), (ii) periodically chirped QPM structures (Bang, Clausen, Christiansen, and Torner 
[1999]), and (iii) optical superlattices (Zhu and Ming [1999]; Fradkin-Kashi, Arie, Urenski, and Rosenman [2002]). 



4.2.1. Phase-reversed QPM structures 

The idea of the phase-reversed QPM structures (Chou, Parameswaran, Fejer, and Brener [1999]) is illustrated in 
Fig. 12(a). Such a structure can be explained as a sequence of many equivalent uniform short QPM sub-structures 
with the length Aph/2 connected in such a way that at the place of the joint two end layers have the same sign 
of the quadratic nonlinearity. Any two neighboring junctions have the opposite signs of the x''^^ nonlincarity. In 
other words, the phase-reversed QPM structure can be thought of as a uniform QPM structure with a change of 
the domain phase by tt characterized by the second (larger) period Aph. Modulation of the quadratic nonlinearity 
d{z) in the phase-reversed QPM structure with the filling factor D = 0.5 can be described by the response fimction 

d{z) = rfo(_l)int(22/AQ)(_-^y„t(2^/Aph)^ (4 6) 



and it can be expanded into the Fourier series, 

+00 +00 +00 +00 



diz) = do E aie'""'^ E ffme'^-^ = do E E 5/™e^'^""^ (4.7) 



/ — — oo m— — oo 1— — 00 m— — oo 



where go = 0, gi^Q ^ i^/nl), g,„^o = (2/7rm), g;™ = gig^, Gi {2tt/Kq)1, = (27r/Aph)m, and dm = 
{2'k/Kq)1 + (27r/Aph)m. 

The phase- reversed QPM structure is characterized by a set of the reciprocal vectors {G/„i}, which are colhnear 
to the normal of the periodic sequence with the magnitude depending on all parameters: I, m, Aq, Aph. Two of these 
vectors can be chosen to phase-match two parametric processes involved into the multistep-cascading, such that 
Giimi = — A5i and Gi^m^ = — A62 ■ Then, the two QPM periods satisfying the double-phase-matching condition 
can be found as follows, 



2TT{l2mi - lim2) 



m2A6i — miA62 



Aph = 



27r(/2mi - /17712) 



;2A6i - ;iA62 



(4.8) 
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Figure 12: Schematic of the non-uniform QPM structures: (a) phase-reversed QPM structure (Chou, Parameswaran, 
Fejer, and Brener [1999]), and (b) periodically-chirped QPM structure (Bang, Clausen, Christiansen, and Torner 
[1999]). 



For the case of the THG multistep cascading process, Eqs. (4.8) are transformed into 

A(^2"^l - hm2) 



(mi — 2m2)n{uj) + 2(toi + m2)n{2uj) — 3min(3a;) 



A. 



ph 



Khmi - lim2) 



(4.9) 



(4.10) 



{h - 2l2)niuj) + 2{li + l2)n{2uj) - 3?i7i(3lj) 

In addition to Eqs. (4.9) and (4.10), the design shown in Fig. 12(a) impose an additional condition that the 
ratio 2Aph/AQ is an integer number. Nevertheless, the corresponding number of the phase-matched wavelengths is 
larger than that achieved in the uniform QPM structure for the collinear geometry. The phase-reverse grating with 
arbitrary ratio 2Aph/AQ was considered theoretically by Johansen, Carrasco, Torner, and Bang [2002] in the study 
of spatial parametric solitons. The study of the effect of arbitrary 2Aph/AQ on the efficiency of different types of 
the frequency conversion processes has not been investigated yet, to the best of our knowledge. Experimentally, 
the phase-reversed QPM structure was employed for the wavelength conversion (Chou, Parameswaran, Fejer, and 
Brener [1999]), and for the cascaded single-crystal THG process (Liu, Zhu, Zhu, Qin, He, Zhang, Wang, Ming, 
Liang, and Xu [2001]; Liu, Du, Liao, Zhu, Zhu, Qin, Wang, He, Zhang, and Ming [2002]). 

4.2.2. Periodically chirped QPM structures 

Periodically chirped QPM structures, presented schematically in Fig. 12(b), have initially been suggested and 
designed for increasing an effective (averaged) third-order nonlinearity in quadratic media (Bang, Clausen, Chris- 
tiansen, and Torner [1999]). We use the terminology "periodically chirped QPM" in the analogy with the chirped 
QPM structures that has a linear growth of the period Aq along the structure. This type of structure is very suitable 
for realizing the multiple phase-matching conditions. The properties of the periodically chirped QPM structures 
have been explored experimentally for achieving a larger bandwidth for the wavelength converters (Asobe, Tadanaga, 
Miyazawa, Nishida, and Suzuki [2003]; Gao, Yang, and Jin [2004]). 
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Figure 13: Basic blocks of the quasi-periodic and aperiodic QPM supcrlattices. 



The periodically chirped QPM structure is characterized by the modulated QPM period Ag that is itself a 
periodic function of z [see Fig. 12(b)]: 

A = Aq +eocos f — j , (4T1) 

Double phase-matching conditions written for this type of structure allow to find two periods of the periodically 
chirped QPM structure; the results are similar to the formulas that describe the phase-reversed QPM structures, 



An = 



Ach = 



27r(/2mi - l\m2) 



7772 A6i — miA62 
27r(/2mi - l\m2) 



(4.12) 



(4.13) 



/2A&1 - ;iA62 

The main advantage of the use of the periodically-chirped QPM structures is the possibility to satisfy the double 
phase-matching conditions for any wavelength in a broad spectral range, i.e. for any pair of A61 and A52. 



4.2.3. Quasi-periodic and aperiodic optical supcrlattices 

Another method of the double phase-matching, studied extensively both theoretically and experimentally, is based 
on the use of the quasi-periodic optical supcrlattices (QPOS) (Zhu and Ming [1999]; Fradkin-Kashi and Arie 
[1999]; Fradkin-Kashi, Arie, Urenski, and Rosenman [2002]) and aperiodic optical supcrlattices (Gu, Zhang, and 
Dong [2000]). In the most of the cases, QPOSs are built with two two-component blocks A and B, as shown 
in Fig. 13, which are aligned in a Fibonacci-like or more general quasi-periodic sequence. Importantly, each of 
the blocks consists of two layers with the opposite sign of the quadratic nonlinearity. In order to illustrate the 
possibility of the double phase-matching in such structures we take, as an example, the structure consisting of two 
blocks aligned in a generalized sequence (Fradkin-Kashi and Arie [1999]). Modulation of the quadratic nonlinearity 
can be described by the following Fourier expansion (Birch, Severin, Wahlstrom, Yamamoto, Radnoczi, Riklund, 
Sundgren, and Wallenberg [1990]): 

ci(z) =c^o^/,„,„e'^™•"^ (4.14) 

where the reciprocal vectors are defined as Gm.n = 27r(?7i -I- nr)/ S, where S = tLa + Lb- 

For phase- matching of two parametric processes with the mismatch parameters A61 and A62, we solve the system 
of equations 

Gm^^m = 27r(mi + niT)/S = -A61, 

(4.15) 

Gm2,n2 = 27r(m2 +n2T)/S = -Afo2, 

and find the corresponding value of 5* and r. The lengths of the blocks La and Lb and a ratio of the sub-layers in 
each block should be found by maximizing fmi.ni and fm2,n2- A resulting designed structure allows simultaneous 
phase- matching in a broad spectral range without constrains on the ratio A62/A61. Equations (4.15) are also valid 
for the Fibonacci-type QPOS lattices but, because the parameter r is fixed, the double phase-matching conditions 
can be satisfied for a limited number of wavelengths (Fradkin-Kashi and Arie [1999]) defined from the equation 
A61/A62 = (toi -I- niT)/(m2 + n2T). 
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Figure 14: Left: 2D nonlinear quadratic photonic crystal composed of a periodic lattice of domains (gray) with the 
reversed sign of x'^"*; and Gb arc the 2D reciprocal vectors. Right: Fabricated 2D hexagonally-poled LiNbOa 
structure with the period 18.05 /im (Broderick, Ross, Offerhaus, Richardson, and Hanna [2000]). 



In several papers Zhao, Gu, Zhou, and Wang [2003]; Gu, Zhang, and Dong [2000], and Gu, Dong, Zhang, 
and Yang [1999] studied aperiodic one-dimensional optical superlattices for the simultaneous phase-matching of 
several parametric processes. In this case, the thickness of the layers and their order can be found by solving an 
inverse problem maximizing the efficiency of both the processes involved into the multistep-cascading parametric 
interaction. Experimental results that employed quasi-periodic and aperiodic optical lattices for realizing the 
multistep parametric interactions are given in Tables 2, 5, and 6. 

4.3. QUADRATIC 2D NONLINEAR PHOTONIC CRYSTALS 

Nonlinear photonic crystals (NPCs) with a homogeneous change of the linear refraction index and a two-dimensional 
(2D) periodic variation of the nonlinear quadratic susceptibility have been suggested by Berger [1998, 1999] as a 
2D analog of the QPM structures. Both theoretical and experimental results published so far show that this kind 
of structures can effectively be employed as a host media for realizing many different types of multistep-cascading 
parametric processes (see, e.g., Broderick, Ross, Offerhaus, Richardson, and Hanna [2000]; Saltiel and Kivshar 
[2000a]; de Sterke, Saltiel, and Kivshar [2001]; Chowdhury, Hagness, and McCaughan [2000]). Below, we discuss 
some of the possible application of these 2D structures for multistep parametric interaction and multi-frequency 
generation. 

4.3.1. Phase-matching in two-dimensional QPM structures 

A schematic structure of 2D NPC is shown in Fig. 14. A simple way to obtain the phase-matching conditions for 2D 
NPC is to use a reciprocal lattice formed by the vectors Ga and Gc, defined as \Ga\ = 27r/Aa and \Gc\ = 27r/Ac . 
For a hexagonal lattice, A^ = Ac = a^/3/2, where a is the distance between the centers of two neighboring inverted 
volumes, the so-called lattice spacing. All reciprocal vectors of the 2D NPC crystal are formed by a simple rule, 
Gm,n = mGc + nGa- Any two vectors of this set can be used to compensate for the bulk mismatch parameters Abi 
and A62, however, the phase-matching conditions require, in most of the cases, non-collinear parametric interactions. 

The diagrams for calculating the phase-matching conditions for the THG multistep-cascading process are pre- 
sented in Fig. 15 for (a) the process uj + tu = 2uj which is phase-matched by the reciprocal vector G,„j_„j , and (b) 
the process uj + 2lo = 3lu, which is phase-matched by the reciprocal vector Gm2,n2- Phase-matching is achieved 
by choosing the lattice spacing a and the angle of incidence f3. 2D NPC can be also used for simultaneous phase- 
matching of three nonlinear processes, e.g. second-, third-, and fourth-harmonic generation (Saltiel and Kivshar 
[2000a]), or the generation of a pair of SH waves and the third- (Karaulanov and Saltiel [2003]) or fourth- (de Sterke, 
Saltiel, and Kivshar [2001]) harmonic generation. Experimentally, the simultaneous second-, third- and fourth- 
harmonic generation was recently observed in 2D poled bulk LiNbOs (Broderick, Ross, Offerhaus, Richardson, and 
Hanna [2000] and Broderick, Bratfalean, Monro, Richardson, and de Sterke [2002]). The first experiment with 
phase-matching in 2D NPC made in a poled LiNbOa waveguide slab was reported in (Gallo, Bratfalean, Peacock, 
Broderick, Gawith, Ming, Smith, and Richardson [2003]). As shown in He, Tang, Qin, Dong, Zhang, Kang, Sun, 
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Figure 15: (a,b) Diagrams of the double phase-matching conditions for single-crystal THG in 2D NPC structures. 
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Figure 16: Multiple second- (0^2), third- {uja), and fourth- (1^4) harmonic generation schemes realized in a 2D NPC 
structure; orange — SHG (766 nm), green — THG (511nm), blue — FHG (383 nm) (Broderick, Ross, Offerhaus, 
Richardson, and Hanna [2000]). 

and Shen [2003], the method of a direct electron-beam lithography can also be used for creating 2D nonlinear 
photonic domain-reversed structures. 

The shapes of the inverted domain is an important property of the 2D NPC structures, and they can be employed 
for an effective optimization of the parametric conversion processes. Optimization of the domain shapes for the 
maximum efficiency of THG and FHG processes was reported by Norton and de Sterke [2003b] . For the case of 
SHG in Lee and Hagness [2003] the following three types of the domain shapes were compared: hexagonal, circular, 
and elliptical (aligned to the direction of SHG). The numerical simulation made by the authors revealed that the 
elliptically-poled domain pattern yields the highest frequency conversion efficiency among these three types of the 
poling structures. The role of the filling factor for the case of SHG was studied both theoretically and experimentally 
in Wang and Gu [2001]; Ni, Ma, Wang, Cheng, and Zhang [2003]. 
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Figure 17: Multi-channel the third- and forth-harmonic generation in 2D NPC structures: (a) the step lo + uj = 2uj\ 
(b) the step ^ + 200 — ^00 for THG; (c,d) the step 2lo + 2w = Alo for FHG; Kp,,, Ki^j, Km,n, and Kfc^; are reciprocal 
wave- vectors of the lattice. 

4.3.2. Multi-channel harmonic generation 

In the first experimental work on the harmonic generation in the 2D NPC structures (Broderick, Ross, Offerhaus, 
Richardson, and Hanna [2000] and Broderick, Bratfalean, Monro, Richardson, and de Sterke [2002]), it was noticed 
that each generated harmonic has multiple outputs at different angles (see Fig. 16). The reason for that is that 
each harmonic can be generated by using several different phase-matching conditions. In several papers, an efficient 
method was suggested to make use of this property and to combine the multiple outputs in order to generate 
efficient harmonics with the additional advantage of being coUinear to the pump. A phase-matching geometry for 
THG in 2DNPC, for which generated TH wave is coUinear to the input wave, was suggested in Karaulanov and 
Saltiel [2003]. This scheme uses the fact that the TH wave is generated though two different channels leading to 
the factor four improvement in comparison with the conventional coUinear scheme. 

A geometry for the multiple phase-matching of this THG multichannel parametric interaction is presented in 
Figs. 17(a,b). The process starts with the generation of a pair of the SH waves by employing the wave vectors k2 
and k2 . The phase-matching conditions are satisfied by two symmetric reciprocal vectors K„i^„ and K„j__„. Each 
SH wave interacts again with the fundamental wave [see Fig. 17(b)] via a pair of the phase-matched interaction with 
participation of the reciprocal vectors Kp ,j and Km._q, thus generating a pair of the coUinear TH waves with one 
and the same wave vector ka . The two TH waves interfere constructively in the direction of the input fundamental 
wave, resulting in the overall higher TH efhcicncy (in the nondepleted regime). 

A scheme for a single-crystal FHG process in a 2D QPM structure where the fundamental and FH waves are 
coUinear was studied in de Sterke, Saltiel, and Kivshar [2001]. The parametric interaction includes cascading of 
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Figure 18: Schematic of the simultaneous optical wavelength interchange in a 2D NPC structure (Chowdhury, 
Staus, Boland, Kucch, and McCaughan [2001]). 

three phase-matched second-order processes that are simultaneously phase-matched, de Sterke, Salticl, and Kivshar 
[2001] demonstrated that the FHG efficiency achieved by this scheme is larger than the comparable schemes by 
a factor that reaches four at low input intensities. The phase- matching geometry for this process is shown in 
Figs. 17(c,d). Two SH waves are generated at the first step [see Fig. 17(a)], and they interact non-coUinearly 
generating a phase-matched FH wave, which itself is coUinear with the fundamental wave [see Fig. 17(c)]. As shown 
in Norton and de Sterke [2003a], there exist two more FHG channels that are simultaneously phase-matched also 
leading to the generated FH waves coUinear to the input fundamental wave [see Fig. 17(d)]. Thus, in this case the 
generated FH wave is produced in result of an constructive interference of three FH waves generated at different 
phase-matching conditions, but propagating in the same direction. 

4.3.3. Wave interchange and signal deflection 

The first experimental demonstration of the simultaneous optical wavelength interchange by the use of a 2D NPC 
structure was reported in Chowdhury, Staus, Boland, Kucch, and McCaughan [2001]. A nonlinear 2D lattice 
fabricated in LiNbOa was designed to provide an interchange of the waves with the wavelengths Ai = 1535nm 
and A2 = 1555nm. The pump was selected at Ap = 777. 2nm. The wavelength interchange process takes place by 
means of the two concurrent difference- frequency generation processes: tUp — uii = L02 and tUp — UJ2 = wi- The two 
difference-frequency processes diffract the converted signals from the unconverted ones, as depicted in Fig. 18. The 
information carried by the input beam 1 will be carried by the output beam 2 and vice versa. In this experiment, 
one of the main characteristics of the 2D NPC phase-matching, i.e. non-coUinearity of the parametric interactions, 
is used as a real advantage in separating two beams at the output. 

Another scheme proposed by Saltiel and Kivshar [2002] is also based on the advantage of non-coUinearity of the 
parametric interactions in the 2D geometry, and it demonstrates how the signal wave can be deflected or split by 
the pump after the interaction. The phase-matching conditions for this scheme arc shown in Fig. 19. The collincar 
pump and signal arc at the same frequency, but they are polarized at the orthogonal directions. If the pump carries 
some information, the signal will be modulated according to this information after the deflection. This interaction 
belongs to the two-color multistep cascading processes discussed above. It can also be considered as a spatial analog 
of the wavelength conversion process discussed above. 

In conclusion of this part of the review, we would like to mention that several theoretical studies predicted that 
the double phase-matched interactions such as the THG multistep parametric process can be realized in photonic 
bandgap (PBG) structures with the efficiency which should be by several orders higher than that of the conventional 
nonlinear media of the similar length. In particular, this includes an infinite PBG system as a host of the second- 
order multistep cascading considered by Konotop and Kuzmiak [1999], cascaded THG in a finite PBG crystal 
discussed in Centini, D'Aguanno, Scalora, Sibilia, Bertolotti, Blocmer, and Bowden [2001], and a photonic-well 
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Figure 19: Phase matching for the first (top) and second (below) cascading steps, yipUip — Z2p, Z2pZis — zis'- 
Drawings correspond to niz < n2z < niy (Salticl and Kivshar [2002]). 

type PBG system analyzed in Shi and Wang [2002]. However, no experimental results on the fabrication of PBG 
structures with strong nonlinear properties were reported yet. 

In conclusion, the part 4 presented a brief overview of different techniques for achieving the simultaneous phase- 
matching of several nonlinear parametric processes in optical structures with a modulated second-order nonlinear 
susceptibility. In all those cases, the double phase-matched interaction becomes possible in a wide region of the 
optical wavelengths provided the QPM structure used to achieve the phase-matching conditions possesses one extra 
parameter (e.g., the modulation period in the chirped QPM structures, or the second dimension, in the case of 2D 
nonlinear photonic crystals). Some of the possible applications of the double phase- matching processes include the 
simultaneous generation of several optical frequencies in a single-crystal structure, multi-port frequency conversion, 
etc. The results presented above look encouraging for experimental feasibility of the predicted effects in the recently 
engineered ID and 2D periodic optical superlattices. 

§ 5. Multi-color parametric solitons 

In the previous sections, we have discussed the features of multiple parametric processes using the plane- wave and 
continuous-wave approximations. However, parametric interactions can strongly modify the dynamics of spatial 
beams or temporal pulses. The parametrically coupled waves in a medium with quadratic nonlincarity may expe- 
rience mutual spatial focusing or temporal compression and lock together into a stationary state, quadratic soliton 
(see Sukhorukov [1988]; Torner [1998]; Kivshar [1997]; Etrich, Ledcrer, Malomed, Peschel, and Peschel [2000]; 
Torruellas, Kivshar, and Stegeman [2001]; Boardman and Sukhorukov [2001]; Buryak, Di Trapani, Skryabin, and 
Trillo [2002], and references therein). 

The first analysis of solitons supported by multistep parametric interactions was reported by Azimov, Sukhorukov, 
and Trukhov [1987] who found that as many as seven waves can be trapped together, and such multi-color solitons 
may be remarkably stable, in particular in the presence of absorption. Since phase mismatch and strength of several 
parametric interactions can be engineered in nonlinear periodic structures (see Sec. 4), there exists a flexibility in 
controlling the soliton properties, making them useful for potential applications including all-optical switching. In 
this section, we overview the properties of multistep parametric solitons that can form under the conditions of 
cascaded THG, two-color FWM, frequency conversion, and FHG. 
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Figure 20: (a,b) Profiles of the THG parametric solitons for different values of phase mismatch (normalized units). 
Labels 1,2, and 3 indicate the number of harmonic. (Kivshar, Alexander, and Salticl [1999]) 



5.1. THIRD-HARMONIC PARAMETRIC SOLITONS 

Generation of the TH wave through cascaded SHG and SPG processes was discussed in Sec. 3.1. Komissarova 
and Sukhorukov [1996] demonstrated that these parametric interactions can result in simultaneous trapping of all 
the interacting waves and formation of a three-colour parametric soliton. Kivshar, Alexander, and Salticl [1999] 
performed a detailed investigation of these solitons and demonstrated their stability under various conditions. Such 
solitons can be described by Eqs. (3.1) with additional terms accounting for beam diffraction. 



dAi ^ i 

dz 2ki 
dA2 i 
^+2fc^""^"^ 
dAs i 

dz 2k^ 



A2 



(5.1) 



where the contribution due to the direct THG process is neglected, assuming that the cascaded x processes are 
dominant. The operator V'^ acts on the spatial dimensions in the transverse plane (perpendicular to the beam 
propagation direction z). Similar equations can describe the formation of temporal solitons, where nonlinearity 
compensates for both the group-velocity mismatch and second-order dispersion effects (Huang [2001]). 

Stationary propagation of solitons is possible only when there is no energy exchange between the constituent 
waves; this requires that the individual phase velocities are synchronized due to nonlinear coupling. Such solutions 
of Eqs. (5.1) have the form A„i(r, z) ~ i3,„(r)e""''^ , where (3 is the nonlinear propagation constant and Bm{r) are 
the transverse soliton profiles. Spatial soliton properties were analyzed in the (l-l-l)-dimensional geometry, when 
the beam is confined by a planar waveguide and experiences diffraction only in one transverse direction so that 
V'^ = d^/dx'^. There exist two families of single-hump solitons, however solutions with higher power are always 
unstable. Soliton examples from the low-power branch arc shown in Pig. 20. It was also found that cascading 
interactions can lead to radiative decay of self-localized beams, however such quasi-solitons can demonstrate robust 
propagation for several diffraction length. 

Lobanov and Sukhorukov [2002, 2003] developed a theoretical description of THG solitons in QPM structures, 
where all the frequency components oscillate along the propagation direction. The spectrum of these oscillations was 
determined analytically and numerically, and it was found that the solitons can be described by averaged equations 
with additional terms accounting for induced Kerr effects: self-phase modulation, cross-phase modulation, and 
third-harmonic generation. 



5.2. TWO-COLOR PARAMETRIC SOLITONS 

As discussed above in Sec. 3.3, two-colour multistep cascading can involve up to six steps corresponding to type I and 
type II interactions between two pairs of orthogonally polarized PP and SH waves. Existence of spatial quadratic 
solitons under these most general conditions was demonstrated by Boardman, Bontemps, and Xie [1998] who found 
that, by changing the PP polarization at the input, the SH output can be precisely controlled. A possibility for 
achieving beam steering and switching based on collision of two solitons was demonstrated as well. 
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Figure 21: (a-c) Characteristic dependences of the total sohton power on the propagation constant corresponding to 
different phase mismatches for two- wave (A-S and B-S) and three- wave (AB-S) sohton famihes (normahzed units). 
Sohd hues show stable solitons, and dashed - unstable; open circles mark the bifurcation points, (the figures 
adapted from Sukhorukov [2000]) 




Figure 22: (a) Change of the normalized power in FF (A, solid) and SH (S, dotted) components, which initially 
constitute a two- wave soliton, and in the guided mode (B, dashed), demonstrating amplification of a guided wave, 
(b-d) Evolution of the intensity profiles for (b) the effective waveguide (SH), and (c,d) A and B FF components, 
respectively. (Kivshar, Sukhorukov, and Saltiel [1999]) 



The general case involving four waves coupled by six parametric processes requires three phase-matching condi- 
tions to be fulfilled simultaneously, and the respective components of the x'^' susceptibility tensor should also be 
nonzero. Formation of three-wave solitons which require two phase-matching conditions may be easier to achieve, 
as indicated by the experimental results in the plane- wave regime (see references in Sec. 3.3). Kivshar, Sukhorukov, 
and Saltiel [1999]; Sukhorukov [2000] considered formation of solitons involving two FF components (A and B) that 
are coupled together through a single SH component (S). The soliton formation can be described by the following 
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set of coupled equations that are obtained in the slowly varying envelope approximation, 




iaBsB e 



<2 iAk-BsZ 



(5.2) 



where (tas and ctbs are proportional to the elements of the second-order susceptibility tensor, as discussed in the 
previous sections, A/cas = — 2fcA and A^bs = — 2/cb are the wave-vector mismatch parameters for the A-S 
and B-S parametric interaction processes, respectively. 

Similar to the case of the THG cascading, multi-component bright solitons are found in the form of stationary 
waves, which phases are locked together. When one of the FF waves is zero, then other two waves can form a 
type I quadratic soliton (A-S or B-S). However, these solitons may experience an instability associated with the 
amplification of the orthogonally polarized FF wave through the parametric decay instability of the SH component. 
Such changes of stability are associated with the bifurcation for two- wave to three- wave solitons, as shown in Fig. 21 
for a (l-l-l)-dimensional case. We note that the Vakhitov-Kolokolov stability criterion (Vakhitov and Kolokolov 
[1973]; Pelinovsky, Buryak, and Kivshar [1995]) can only be used as a necessary condition: branches with the 
negative power slope are unstable, but the positive slope does not guarantee stability. 

Examples in Fig. 21 demonstrate the existence of stable solitons that have the same power but different polar- 
izations of the FF waves. Such multistability is a sought-after property on nonlinear systems, since this may allow 
realization of controlled switching between different states. In Fig. 22, we illustrate the development of the soliton 
instability which results in a power exchange between the two FF components. Such type of polarization switching 
was earlier predicted for plane waves by Assanto, Torelli, and Trillo [1994]. 

The properties of two-color multistep cascading solitons were analyzed by Towers, Sammut, Buryak, and Malomed 
[1999]; Towers, Buryak, Sammut, and Malomed [2000] for the case when two FF waves and the SH component are 
additionally coupled together through the type II parametric process. Due to this interaction, all solitons contain 
both FF waves, and transition between the states with different polarizations along the soliton family does not 
involve bifurcations. It was found that soliton multistability can be realized under appropriate conditions, both in 
planar (one-dimensional) and bulk (two-dimensional) configurations. On the other hand, Towers, Buryak, Sammut, 
and Malomed [2000] reported that no multistability occurs if the BB-S interaction is suppressed and the waves are 
coupled only through A A-S and AB-S processes. 

5.3. SOLITONS DUE TO WAVELENGTH CONVERSION 

Formation of spatial parametric solitons requires a strong coupling between the interacting waves, and this can be 
achieved when the parametric interactions are nearly phase-matched. In the cases of parametric THG and two- 
colour FWM, two (or even more) phase-matching conditions should be satisfied simultaneously. On the other hand, 
in the wavelength conversion scheme (see Sec. 3.2), only single phase-matching condition should be implemented, 
and this greatly simplifies the requirements for experimental observation of spatial beam localization and soliton 
formation under the conditions of multistep cascading. Indeed, the first experimental study was recently reported by 
Couderc, Lago, Barthelemy, De Angelis, and Gringoli [2002] who demonstrated that a multi-color soliton composed 
of a pump and its second harmonic creates an effective waveguide that can trap a weak probe beam which frequency 
is slightly detuned from the pump wave (see Fig. 23) . This trapping is realized due to a parametric coupling of the 
probe with both the pump and harmonic waves, and additional side-band frequencies are generated in the process 
according to the principles of wavelength conversion. Couderc, Lago, Barthelemy, De Angelis, and Gringoli [2002] 
demonstrated that evolution of a weak probe is governed by the equation. 



where A and B arc the profiles of mutually trapped pump wave and its second harmonic components, 6uj is the 
frequency detuning of the probe, loi and ki are the frequency and the wavenumber of the fundamental frequency 
wave, and Afc is the phase mismatch between the SH and FF components. The last term in Eq. (5.3) defines the 
profile of an effective waveguide which is experienced by the probe beam, and trapping can be realized when the 
overall sign is positive. It has been demonstrated in earlier studies that the FF component (|A|^) dominates in a 
quadratic soliton when Afc > 0, whereas the SH component (|-B|^) becomes larger at negative mismatches, Afc < 0. 




(5.3) 
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Figure 23: Beam profiles at the output of the 2cm long KTP crystal for low and high input intensities. Horizontal 
scale: wavelength; vertical scale: transverse dimension on the output pattern. Right part of the figure corresponds 
to the probe wave (1061 nm), the left part corresponds to the FF component of the soliton (pump at 1064 nm). 
(Couderc, Lago, Barthelemy, De Angelis, and Gringoli [2002]) 
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Figure 24: (a) Thick — two- wave (SH + FH), and black — three- wave solitons; solid and dashed lines mark stable 
and unstable solutions, respectively. Open circle is the bifurcation point; (b) Development of a decay instability of 
a two- wave soliton corresponding to /3 = 1 in plot (a), and generation of a three-component soliton. Dotted, solid, 
and dashed curves in the left plot show the FF, SH, and FH peak normalized intensities vs. distance, respectively. 
(Sukhorukov, Alexander, Kivshar, and Saltiel [2001]) 

Therefore, the probe can be trapped in both cases, however anti-wavcguiding effect may occur for particular values 
of the mismatches and pump intensities. 

5.4. OTHER TYPES OF MULTI-COLOR PARAMETRIC SOLITONS 



Formation of spatial parametric solitons due to FHG in a planar waveguide configuration was analyzed theoretically 
by Sukhorukov, Alexander, Kivshar, and Saltiel [2001] . It was found that, similar to the case of two-color parametric 
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solitons (see Sec. 5.2), there can exist two-wave solitons (coupled second and fourth harmonics) and three-wave 
soUtons containing aU the three frequency components. Such coexistence of different solitons may give rise to 
multistability, as illustrated in Fig. 24(a). A two- wave soliton can exhibit parametric decay instability due to the 
energy transfer into the fundamental-frequency wave, this process is illustrated in Fig. 24(b). 

Towers and Malomcd [2002] performed a theoretical analysis of polychromatic solitons consisting of two low- 
frequency components uJi and uj2 and three high-frequency waves, 2wi, 2w2; and uji + uj2- It was found that the 
power of a polychromatic soliton can be smaller compared to conventional two- wave parametric solitons, and this 
can be advantageous for applications. It was also demonstrated that polychromatic solitons can emerge after a 
collision of two-wave solitons, and soliton interactions may be used to implement a simple all-optical XOR logic 
gate. 

§ 6. Conclusions 

Parametric interactions and phase matching are the key concepts in nonlinear optics. The generation of new 
waves at the frequencies which are not accessible by standard sources and the efficient frequency conversion and 
manipulation are among the main goals of the current research in nonlinear optics. The QPM technique is becoming 
one of the leading technologies for the optical devices based on the parametric wave interaction; it allows generating 
new harmonics and can be made compatible with the operational wavelengths of optical communication systems. 

In this paper, we have presented, for the first time to our knowledge, a systematic overview of the basic principles 
of the simultaneous phase matching of two (or more) parametric processes in different types of one- and two- 
dimensional nonlinear quadratic optical lattices, the so-called multistep parametric interactions. In particular, we 
have discussed different types of multiple phase-matched processes in the engineered QPM structures and two- 
dimensional nonlinear quadratic photonic crystals, as well as the properties of multi-color optical solitons generated 
by the multistep parametric processes. We have also summarized the most important experimental demonstrations 
for the multi-frequency generation due to multistep parametric processes. We believe that such a comprehensive 
summary of the up-to-date achievements will become a driving force for the future even more active research in 
this exciting field of nonlinear optics. 
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